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INHOMOGENEOUS NAVIER-STOKES EQUATIONS IN THE 
HALF-SPACE, WITH ONLY BOUNDED DENSITY 

RAPHAEL DANCHIN AND PING ZHANG 

Abstract. In this paper, we establish the global existence of small solutions to the inhomo- 
geneous Navier-Stokes system in the half-space. The initial density only has to be bounded 
and close enough to a positive constant, and the initial velocity belongs to some critical 
Besov space. With a little bit more regularity for the initial velocity, those solutions are 
proved to be unique. In the last section of the paper, our results are partially extended to 
£SJ , the bounded domain case. 
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Oh ' 1- Introduction 

We are concerned with the global well-posedness issue for the initial boundary value prob- 
lem pertaining to the following incompressible inhomogeneous Navier-Stokes equations: 

dtp + div (pu) = 

dt(pu) + div(pu <g) u) — fJ.Au + VII = 

(1.1) I div u = 
u = 

^ ! [ p\t=o = Po, pu\ t =o = pou 

where p = p(t, x) € M. + , u = u(t, i)el and II = II(t, i)eR stand for the density, velocity 

field and pressure of the fluid, respectively, depending on the time variable t € M+ and on 

_il ■ the space variables x £ Q. The positive real number p stands for the viscosity coefficient. 

We mainly consider the case where Q is the half-space M + , except in the last section of the 

paper where it stands for a smooth bounded domain of M (d > 2). 

The above system describes a fluid that is incompressible but has nonconstant density. 

Basic examples are mixture of incompressible and non reactant flows, flows with complex 

structure (e.g. blood flow or model of rivers), fluids containing a melted substance, etc. 
H ■ 

A number of recent works have been dedicated to the mathematical study of the above 

system. Global weak solutions with finite energy have been constructed by J. Simon in [19] 

(see also the book by P.-L. Lions [17] for the variable viscosity case). In the case of smooth 

data with no vacuum, the existence of strong unique solutions goes back to the work of O. 

Ladyzhenskaya and V. Solonnikov in [15]. More recently, the first author [8] established the 

well-posedness of the above system in the whole space M. in the so-called critical functional 

framework for small perturbations of some positive constant density. The basic idea is to use 

functional spaces that have the same scaling invariance as (1.1), namely 

(1.2) (p, u, U)(t, x) i — > (p, Xu, A 2 II)(A 2 t, Ax), (p , u )(x) i — > (p , Xu )(Xx). 
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2 R. DANCHIN AND P. ZHANG 

More precisely, in [8], global well-posedness was established assuming that 

||/9o — 111 d + p~ l ||uo || d <1. 

B^ x (R d )nL°°(R d ) B^i(R d ) 

Above Bp r (R ) stands for a homogeneous Besov space on R (see Definition 2.1 below). 
This result was extended to more general Besov spaces by H. Abidi in [1], and H. Abidi and 
M. Paicu in [2], and to the half-space setting in [10]. The smallness assumption on the initial 
density was removed recently in [3, 4]. 

. d 

Given that in all those works the density has to be at least in the Besov space Bp )OC , one 

cannot capture discontinuities across an hypersurface. In effect, the Besov regularity of the 

i 

characteristic function of a smooth domain is only Bp )!X) . Therefore, those results do not 
apply to a mixture of two fluids with different densities. 

Very recently, the first author and P. Mucha [11] noticed that it was possible to establish 
existence and uniqueness of a solution in the case of a small discontinuity, in a critical 
functional framework. More precisely, the global existence and uniqueness was established 
for any data (pq,uq) such that for some p £ [1, 2c?) and small enough constant c, we have 

(1-3) l|Po-l|| - 1+ d +/i~ 1 ||u || _ 1+ d < c. 

M{B vl p (R d )) B pl ?(R d ) 

— 1+- 
Above, || • || _i + d is the multiplier norm associated to the Besov space B 1 v (R ), 

M(B pl p (R d )) ' P , 

which turns out to be finite for characteristic functions of C 1 domains whenever p > d — 
1. Therefore, initial densities with a discontinuity across an interface may be considered 
(although the jump has to be small owing to (1.3)). As observed later on in [12], large 
discontinuities may be considered if the initial velocity is smoother. In fact, therein, any 
initial density bounded and bounded away from is admissible. Let us emphasize that 
in both works ([11] and [12]), using Lagrangian coordinates was the key to the proof of 
uniqueness. 

A natural question is whether it is still possible to get existence and uniqueness in a critical 

functional framework where po is only bounded and bounded away from zero. As regards 

existence, a positive answer has been given recently by J. Huang, M. Paicu and the second 

author in [14], in the whole space setting, and uniqueness was obtained if assuming slightly 

more regularity for the velocity field. Let us emphasize that once again using Lagrangian 

coordinates is the key to uniqueness. Therefore, assumptions on the initial velocity have 

to ensure the velocity u to have gradient in Lj- (R+ ; L°° (M )) in order that Eulerian and 

Lagrangian formulations of the system are equivalent. While this property of the velocity 

— i+— ■ — i+- 

field holds true if uq is in B l v (R ) , it fails if uq is only in B p<r p (R ) for some r > 1 . As 

a matter of fact, the question of uniqueness in a critical Besov framework for the velocity is 

open unless r = 1 (this latter case requires stronger assumptions on the density, as pointed 

out in [11]). 

In the present work, we aim at extending the results of [14] to the half-space setting. 
Because we shall consider only perturbations of the reference density 1, it is natural to set 



(1.4) 
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a = 1/p — 1 so that System (1.1) translates into 

d t a + u • Va = in R + xft, 

d t u + u- Vn + (l + a)(Vn-/xAn) = in M+ xfi, 

div n = in M + xft, 

n = on R + x<9ft, 

(a,u)\ t =o = (oo,u ) in ft. 

As in the whole space case considered in [14], the functional framework for solving (1.4) is 
motivated by classical maximal regularity estimates for the evolutionary Stokes system. In 
effect, the velocity field may be seen as the solution to the following system: 

(1.5) d t u-nAu + VII = -u- Vu + a(/j,Au- VII), divn = 0. 

In the whole space case, we have for any 1 < p, r < oo and t > 0, 

||(Stu,/iVVVn)|| L r (LP ) < C(/x 1- r||n || .2-2 + \\u-Vu\\ L r {LP) + ||a(/iAu- Vn)|| L r (I , P ) 



def 



~ 



L r ((0,t);LP(R d ))- 



where we have used the notation ||2||£,»-(lj>) 

Given that ||a(£)||x,oo = ||ao ||z,°° f° r all time, it is clear that the last term may be absorbed 

by the l.h.s. if ||ao||z,°° is small enough. Now, it is standard (see Corollary A.l below) that 

in the simpler case where n just solves the heat equation with initial data no, having An in 

■ — 1+- 
L r (0,T; LP) is equivalent to no €. B P)T p provided —1 + ^ = 2 — ^- This implies that p and 

r have to be interrelated through p = 3^2, and thus 3 < p < d. 

Here we aim at extending this simple idea to the half-space setting, or to C 2 bounded 
domains. 

In the half-space case, according to the above heuristics and because homogeneous Dirichlet 
boundary conditions are prescribed for the velocity, the natural solution space for (n, VII) is 

X% rd = {(n,Vn) with nGC([0,T];,Bp,; f (]R^)) and d t u, V 2 n, VII € L r {0,T; L p + )\, 



2 



where L p + = L P (R + ) denotes the Lebesgue space over M + , and B P: r r (M + ) stands for the set 



2 



of divergence free vector fields on R + with Besov regularity B p<r r (M + ) and null trace at the 
boundary (the exact meaning will be given in Definition 2.3 below). 

We also introduce the following norm for all T > 0: 

(1.6) ||(n, VII)|| xr =^\H\ LOO{OT . K -^ K)) + ||(^n, M V 2 n,Vn)|| L , (0iT;LP+) , 

and agree that X p ' r and || • ||xp> r correspond to the above definition with T = +00. 

Before stating our main results, let us clarify what we mean by a weak solution to (1.4): 

Definition 1.1. A global weak solution of (1.4) is any couple (a,u) satisfying: 
• for any test function <p € C£°([0, 00) x ft), there holds 



oc 



a(dt<f> + n • V0) dxdt + / <f>(0, -)ao dx = 0, 
'-, -\ Jo Jn Jn 

/ cf> div n dxdt = 0, 
Jn 
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' for any vector valued function <J> = ((J* 1 , • • • , <3? d ) G C£°([0, oo) x Q) , one has 

/ / \u-d t <5> + u®u : V$ + (l + a)(uAu- Vn)$)dxdt + / u ■$((), -)dx = 0. 
Jo isi^ J ic 



Our main statement reads: 



-1+- , ■ — n d 



Theorem 1.1. Let a € L~ and u € 6 p , r p (R*) D B- r p (R|) with p= ^ <p < J^ 

and r G (1, oo). There exist two positive constants c$ = co(r, d) and c\ = ci(r, d) so that if 

(1.9) 1 1 i*o [ | .-i+d < c /i and ||a ||L°f < ci 

B„ r (IRj.) 



r/ 



then (1.4) has a global solution (a,u, VII) in the meaning of Dehnition 1.1 with f2 — ^ + , 
satisfying ||a(£)||i°° = ||ao |[z,^° for all t > and u G X p ' r n I p,r . Moreover, there exist 
Ci = Ci (r, d) and C2 = C2 (r, d) so that 



(1.10) 


ll(« ; 


,vn)|| X P 


i__L 

,r + / Lt ± 2r 


and, if 


j_ def 1 _ 


r-l 




(1.11) 


IKu.vnjiixp.r + M 1 "^ 


If in addition p > d then 


Vu € L 9 C 


(1.12) 









|Vu|| dr < Clfl r \\uo 



-1+ 



d 



L^(R +; lF^) B v ,r P (Rl) 



— 1 — - 

2r IIVtl|| L 2r( R+;L a) < C2fi r \\U0\\ , , ll 



B 



H 



fi 



; L^°) with g = d+{ l p _ 1) ~ and 



L9(R,;L°°) < C2 1 1 1*0 1 1 ._l+d 1 

B-, + p (R d + ) 



and uniqueness holds true. 

Remark 1.1. In contrast with the whole space case in [14, 18], it is not clear that one may 
improve the isotropic smallness condition (1.9) to an anisotropic one allowing for arbitrarily 
large vertical velocity. The reason why is that even for solutions to the linear Stokes system 
in the half-space, the horizontal component of the velocity u h depends on both Uq and Uq 
(see the formula in Theorem 2.1 below). 

Remark 1.2. We shall extend this statement to a more general critical (or almost critical) 
Besov setting, see Theorems 5.1 and 5.2 below. We could also deal with the local well- 
posedness of (1.4) with general velocity and small inhomogeneity. For a clear presentation, 
we skip the details here. 

Let us briefly describe the plan of the rest of the paper. The next section is devoted to the 
linearized velocity equation of (1.4) in the half-space, that is the evolutionary Stokes system. 
We first derive an explicit solution formula in the spirit of that of S. Ukai in [20], and then 
deduce maximal regularity type estimates similar to those of the whole space. We consider 
the general situation with prescribed (possibly nonzero) value for divu as it will be needed 
when reformulating (1.4) in Lagrangian coordinates. The next two sections are devoted to 
the proof of the existence part of Theorem 1.1, first under a stronger assumption on the 
density, and next in the rough case corresponding to the hypotheses of the theorem. The 
case of more general Besov spaces will be examined in Section 5. The proof of uniqueness 
is postponed in Section 6. In the final section, we partially generalize Theorem 1.1 to the 
bounded domain setting. Some technical lemmas related to maximal regularity and L p (L q ) 
estimates for the heat equation in the whole space (or Stokes system in bounded domains) 
are presented in Appendix. 
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2. The evolutionary Stokes system in the half-space 
This section is devoted to the study of following system in the half-space: 



(2.1) 



' d t u - fiAu + Vn = / in R+ x R^, 

divu = g in K + xKf, 

n = on R+x<9R^, 

u\t=o = ^0 on ^+ • 



We shall first derive an explicit formula for the solution to this system, and next prove the 
key a priori estimates that are needed for getting the main results of our paper. 

2.1. A solution formula. This part extends a prior work by S. Ukai [20] (see also [7]) to 
the case where there is a source term / in the velocity equation, and where the divergence 
constraint is nonhomogeneous. Let us recall that in [20], it was assumed that / = and 
g = (but u need not be zero at the boundary), and that in [7] nonzero / was considered 
(but still u is divergence free and / has trace zero). Furthermore, the gradient of the pressure 
was not computed therein, a computation that turns out to be essential for us as VII appears 
in the right-hand side of the velocity equation (1.5). 

Before writing out the formula, let us introduce a few notations. We denote A = Yle=i ®xe 
and Ah = X^=i ^L> anc ^ define \D\ ±l and l-D/il^ 1 to be the Fourier multipliers with symbol 






l&l = >At and l£l = \yji) > respectively. 



The notations Rj and Sj stand for the Riesz transforms over R and R fc ~ , namely 
Rj = dj\D\~ l for j = !,••• ,d and Sj = dj\Dh\~ l for j = 1, ■ ■ ■ ,d — 1. 

We further set R^ = (R±,--- , i^-i) and S = (Si,--- ,S d -i). 

As in [20], for u = (u h ,u d ) with u h = (u 1 , - - • , u d_1 ), we define the operators V d and Vh 
by 

(2.2) V d u = f -S ■ u h + u d and V h u A =u h + Su d . 

We shall see later on that both V^u and V d u satisfy a heat equation, this is the main 
motivation for considering those two quantities. 



We also denote x = (xh,x d ) with Xh = (x±, • • • ,x d -i)- Let r be the restriction operator 
from R to R + , that is rf = /| K d , and eo(/),e a (/),e s (/) be the extension operators given 
by 

p ( f \ - / / for x d > °; . (f \ _ / f(x) for x d > 0, 

° U; 1 for x d <0, Ba{1) \ -f(x h ,-x d ) for x d < 0, 
(2.3) 



i /• r\ f f (#) for x d > 0, 

and e s (f) = < J ,) ' x , ' 

V ' \ f(x h ,-x d ) for x d <0. 



When solving (2.1), we shall repeatedly consider the following two equations: 

(2.4) (d d +\D h \)w = \D h \f in M.%, and -fw = on dR^_ . 

(2.5) (d d + \D h \)v = f in R^_, and 7 t> = on <9R^. 
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Here and in what follows, 7 stands for the trace operator on <9K + . 

Finally we denote by H the harmonic extension operator from the hyperplane dR + to 
the half-space M + . More precisely, for any given b : <9M + — > R, we set Hb to be the unique 
solution going to zero at 00 of 

(2.6) AHb = in R*, and ^Hb = b on SR*. 

Introducing the Fourier transform Th with respect to the horizontal component Xh, the 
function ^(Hb) is explicitly given by the formula 

(2.7) Fh(Hb)(Zh, Xd) = e-^ Xd T h b(^ h ), £ h G R d ~\ x d > 0, 
hence in particular 

(2.8) (d d +\D h \)Hb = in R* . 

Lemma 2.1. For any smooth enough data f decaying at 00, Equation (2.4) has a unique 
solution going to at 00, and Equation (2.5) has a unique solution with gradient going to 
at infinity Furthermore, denoting by U and P the solution operators for (2.4) and (2.5), 
one has 



Vd) dy d 



(2.9) Uf = rR h -S{R h -Se a {f) + R d e s (f)) and Pf(x h ,x d ) = Al - U)f{x h , 

Jo 

and the following identities are satisfied: 

(1) V h U = UV h ; 

(2) d d U = (I-U)\D h \: 

(3) V h P = SU; 

(4) 8 d P = I-U; 

(5) AP = d d -\D h \; 

(6) [P,d d ] = -H 7 . 

Proof. If w is a solution to (2.4) then it also satisfies the following Poisson equation: 

-Aw = (\D h \-d d )\D h \f 
w\ Xd=0 = 0, 

the unique solution (decaying to at infinity) of which is given by 

w = r(-A)- 1 e a ((\D h \-d d )\D h \f). 

As e a (\D h \ 2 f) = \D h \ 2 e a (f) and e a (d d \D h \f) = d d (e s (\D h \f)) = d d \D h \e s (f), we get the 
formula for Uf. 

It is obvious that U commutes with V/j. As regards commutation with d d , we notice that, 
by definition of Uf, 

d d Uf+\D h \Uf=\D h \f, 
which yields (2). 

It is clear that \D h \Pf satisfies (2.4), hence \D h \Pf = Uf and V h Pf = SUf. Similarly, 
the equation for Pf yields 

d d Pf = f-\D h \Pf = f-Uf, 

hence integrating with respect to the vertical variable gives the expression for P. 
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The next item is a direct consequence of the definition of P (just apply dd — \D^\ to the 
equation). Finally, we have by definition of Pddf, 

(d d + \D h \)Pd d f = d d f = (d d + \D h \)f - \D h \f. 

Hence, using (2.8), 

(d d + \D h \)(Pd d f-f + Hjf) = -\D h \f and (Pd d f -/ + H>yf)\ Xd=0 = 0. 

Therefore, using the definition of U, one may write 

Pd d f -f + H 1 f = -Uf. 

Because d d P = I — U, it is easy to complete the proof of the last item. □ 

Remark 2.1. For functions vanishing at x d = 0, operator U coincides with the expression 

rR h -S(R h -S + R d )e 

that has been introduced in [20] and plays the role of the left-inverse of (Id-\-\Di l \~ 1 d d ) for 
functions of K + vanishing at <9R + . Our definition of operator U is slightly more general as 
it allows us to consider functions that do not vanish at x d = 0. 

The main result of this subsection reads: 

Theorem 2.1. Given smooth and decaying data uq, f and g with g = divQ, the unique 
solution (u, VII) of (2.1) is given by 



u =re 



fitA ( 



Whuo) - fiSPg - SUe» tA e a {V d u Q ) 
ft ft 



( 2 10) - SU I e^- T ^ A e a (Nf + Gk)dr + r I e^- T ^e a {Mf + SGk) dr, 

u d =nPg + Ue» tA e a (V d u ) + U I e^- T ^ A e a (Nf + Gk) dr, 

Jo 



and 



(2.11 



V h II = r(R h - S(SR d + R h )) (S ■ e a (f h ) + e s (f d )) + S(U - I)Gk 

+ fi(V h - Sd d )g + Sd t {S ■ UQ h + (/ - U)Q d - H 7 Q d ) + SUNf 

+ (r(\D h \ - 8 d )V h + SUA) (e» tA e a (V d u ) + f e^^ A e a (Nf + Gk) dr) , 

d d U = f d + fi(d d - \D h \)g -d t (S- UQ h + (I - U)Q d - H 7 Q d ) - U(Nf + Gk) 
- (r(|£>„| - d d )\D h \ + UA) (e» tA e a (V d u ) + j e^-^ A e a {Nf + Gk) dr 

where Gk, Mf and Nf are given by 

Gk = -r{R d -R h -S){R h -e a {k h ) + R d e s {k d )) with k = d t Q-fiVg, 

(2.12) Nf = r{ [1 + R 2 d - R d R h ■ S}e s (f d ) + R 2 d S ■ e a {f h ) + R d R h ■ e a (f h )}, 

Mf = rS[R d -R h -S] (R h ■ e a (f h ) + R d e s (f d )) + V h f. 
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Proof. We shall essentially follow the arguments in [7, 20]. Note that setting 

(2.13) u new (t,x) = iMu oU (irH,x), U new (t,x) = n o i d (n~H,x), f new {t,x) = foidi^^x) 

reduces the study to the case fj, = 1, an assumption that we are going to make in the rest of 
the proof. 

The basic idea is to reduce the study to that of the heat equation for the auxiliary functions 
VhU and V d u. As a first step, let us compute II in terms of div /, g and of its trace at 9R + . 
Taking space divergence to (2.1) yields 



-An = d t g -Ag- div / in ™ d 



7 n = 6 on dP d 



II — )■ as |x| — > oo, 



the solution of which is given by 

U = r(-A)- 1 e a (div(k-f)) + Hb with k = d t Q-Vg, 

which along with (2.8) implies that 

, 01/1 , (d d + \D h \)U = - r (^ + |^|)(-A)- 1 e a (div(/-A ; )) 

(2.14) 

= Mea(div(/-fc)) 

with 

(2.15) Mh A M-r(d d +\D h \)(-A)- l h. 

Let z = f {d d + \D h \)u d and Nf d = (d d + \D h \)f d - d d Me a (dw /). We infer from (2.1) and 
(2.14) that 

d t z -Az = Nf + d d Me a (div k) in M + xRf, 
(z-g)\ Xd=0 = 0, 

z\t=o = (d d +\D h \)4. 

Note from the definition of M in (2.15) that 

(2.16) d d Me a (h) = e a (h) - r\D h \(d d +\D h \)(-A)- l [e a (h)]. 
Therefore, because 

(2.17) e a (div k) = div h e a (k h ) + d d e s {k d ) 
and 

(2.18) Me a (divk) = Gk, 
we get 

Nf=\D h \f d -rdiy h e a {f h ) + r\D h \{d d + \D h \)(-A)- l [d\w h e a {f h ) + d d e s {f d )] 

= \D h \Nf. 

Now, using (2.16) with h = dtg — Ag = div A; and (2.17), we thus obtain 

d t (z-g)-A(z-g) = \D h \(Nf + Gk) in l + xlf, 
(z-g)\ Xd =o = in IR+ x9M^, 

(z - g)\t=o = \D h \V d u in R d + . 
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Taking advantage of the solution formula for the heat equation in M + with homogeneous 
Dirichlet boundary conditions, we deduce that 

(z - g)(t) = r\D h \(e tA e a (V d u ) + J e^ A e a (Nf + Gk) dr 

As z — g = (l-D/il + d d )(u — Pg) and u — Pg vanishes at x d = 0, keeping in mind the 
definition of U, we get the second equality of (2.10). 

To derive the solution formula for u , we look at the equation satisfied by VhU. Thanks 
to (2.1), (2.12) and (2.14), we get, observing 

Vhf ~ S(\D h \ + d d )U = V h f - SMe a {d\w{f - k)) = Mf + SGk, 

that 

d t V h u - AV h u = Mf + SGk in M + xR{, 
V h u\ Xrl=0 = in M + xR d + , 



V h u\ t =o = V h u in id 



so that 



(2.19) V h u = re tA e a (V h u ) + r f e^-^ A e a {Mf + SGk) dr. 

Jo 

Because u h = V^u — Su d , combining the above identity and the second formula of (2.10), we 
obtain the solution formula for u h . 

Let us finally derive (2.11). By virtue of (2.1) and (2.10), we may write 

d d U = f d - (dt - A)u d 

= f d - (8 t - A)Pg - (8 t - A)u(e tA e a (V d u ) + J e^ A e a (Nf + Gk) dr) . 

However, because 

AUh = r(d d -\D h \)\D h \h, 

one may write 

(dt - A)Uh = U(d t - A)h + r(\D h \ - d d )\D h \h + UAh. 

Note also that, by virtue of Lemma 2.1, 

d t Pg = d t (Pdw h Q h + Pd d Q d ) = d t (S- UQ h + (I-U)Q d - HjQ d ) 

and that 

APg = (d d -\D h \)g, 

which ensures that 

(2.20) d d U = f d - d t (S ■ UQ h + (I- U)Q d - H 7 Q d ) + (d d - \D h \)g - U(Nf + Gk) 

- (r(\D h \ - d d )\D h \ + UA) (e tA e a (V d u ) + J e^ A e a (Nf + Gk)) dr) ■ 

On the other hand, by virtue of (2.14), one has 

\D h \U = -d d U + Me a (div(/ - k)), 
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and 

Me a (div /) = -(d d + IAJX-A)- 1 (d\v h e a (f h ) + d d e s (f d )) 

= e s {f d ) - (d d + \D h \)(-A)- 1 {dw h e a (f h ) + \D h \e s (f d )), 
which together with (2.20) gives rise to (2.11). This completes the proof of Theorem 2.1. □ 

2.2. A priori estimates. Let us first briefly recall the definition of homogeneous Besov 
spaces in R . Let \ '■ ^ — > [0, 1] be a smooth nonincreasing radial function supported in 
B(0, 1) and such that x = 1 on B(0, 1/2), and let 

The homogeneous Littlewood-Paley decomposition of any tempered distribution uonl is 
defined by 

A k u d =<p(2- k D)u = ^- 1 (<^(2- fe -)^), keZ 

where J- stands for the Fourier transform on R . 

Definition 2.1. For any s£l and (p,r) 6 [l,+oo] 2 , the homogeneous Besov space -B*(R ) 
stands for the set of tempered distributions f so that 

Wf\\B^ r (R d ) = lr H A fc/llLP(R d )llr-(z) < °° 
and for all smooth compactly supported function 9 over R , we have 

(2.21) lim 0(\D)f = in L°°(R d ). 

Remark 2.2. Condition (2.21) means that functions in homogeneous Besov spaces are re- 
quired to have some decay at infinity (see [5] for more details). In particular, we have 

(2.22) / = £>*/ in S'(R d ) 

kez 
whenever f satisfies (2.21). In this paper, we will only consider exponents s < d/p so that 
for f with finite Bp r (R ) semi-norm, (2.21) and (2.22) are equivalent. 

The homogeneous Besov spaces on the half-space are defined by restriction: 

Definition 2.2. For any s G R, and (p,r) € [1, +oo] 2 we denote by B*(R^.) the set of 
distributions u on R + admitting some extension -u € i?f r (R ) on R . Then we set 

ii ii '^ ■ f 1 1 ~ ii 

H U H% ir (R^) = mt W U \\B^ r (R d ) 

where the infimum is taken on all the extensions of u in Bp r (R ). 

We also need to introduce some spaces of divergence free vector fields vanishing at the 
boundary d R + . We proceed as follows: 

Definition 2.3. For 1 < p < oo and < s < 2, we denote by ,6p r (R + ) the completion 
of the set of divergence free vector fields with coefficients in W 2 ' P (R. + ) n Wg' p (R + ) (where 
W 1,P (R+) stands for the subspace of V4 /1,P (R+) functions with null trace at <9R+) for the 
norm || • ||^ r (R^)- 

It is classical (see e.g. [10]) that spaces (Bt, r (R + )) rf (with the divergence free condition) 
and Z3p r (R + ) coincide whenever 1 < p, r < oo and < s < 1/p. 
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The following result extends Lemma 3.2 of [20] to the context of Besov spaces. 

Lemma 2.2. Operators R^, R^ and S map L P (M. ) in itself for any 1 < p < oo, and, with 
no restriction on s,p,r, we have 

||^z||^ r(Rd) < C||z||^ (Rd) for K G {R h , R d , S}. 

Operators Vh, Vd, U, G, M and N map L p + in itself if 1 < p < oo, and B^ r (R + ) in itself 
if 1 < p, r < oo and < s < 2. 

Proof. The result in the Lebesgue spaces just follows from the fact that all those operators are 
combinations of Riesz transforms so that Calderon-Zygmund theorem applies. The result in 
homogeneous Besov spaces stems from the fact that the Riesz operators are Fourier multipliers 
of degree 0, hence map any homogeneous Besov space in itself. □ 

We are now ready to establish a hrst family of a priori estimates for System (2.1). 

Proposition 2.1. Let 1 < p, r < oo and the data uq, f, g fulfill the following hypotheses: 



p,'. 

p 



n € Bp,r r (Ml 



. feL r (R + ;L p +/ 

• V<? € L r (lR + ; L p + ), g = divQ with dtQ € L r (R + ; L p + ) and the following compatibility 
conditions are fulfilled 1 

(2.23) 7^0 = 0, <7|t=o = 0, and 6\( 7 Q d ) = 0. 
Then System (2.1) has a unique solution (u, VII) with 

MeC 6 (l + ;Bj,7(4)) and d t u,V 2 u,VU£L r (R + ;L p + ), 
with also Vu £ L 9 (IR + ;L™) whenever q £ [r, oo) and m € [p, oo] satisfy 

2 2 d , d d 2 2 

0<1 h-<- and — = --H 

r q p m p r q 

Furthermore, the following inequality is fulfilled for all t > : 

(2.24) w ^|M| L ^_a (Rd+)) +^ + "||Vn|| L?(L!F) + \\(d t u,pV\VU)\\ L r {LP+) 



< C(V^|K|| 2 2 + \\(f,n.Vg,d t Q)\\ L r {LP \ 



Bp,r T (R+) 



^2-- 



Finally, if g = then we have u € C(M+; B v y (M+)). 

Remark 2.3. As regards the bounds for Vu, we shah often use the following two cases: 

• P= 5F=2, q = 2r andm= ^, 

• p > d, q = - ■ , g 2 - . and m = +oo. 

Proof. We concentrate on the proof of the estimates in F r (0, T; L^_) for V 2 u and VII. Indeed, 
once the pressure has been determined, u may be seen as a solution of the heat equation 
with source term in L r (R + ; L p + ), the solution of which is given by 

(2.25) u(t) = r (V A e a (n ) + f e^- T ^ A e a (f - VII) dr) ■ 

■2— - 

The condition on uq is ensured by the fact that «o € B v , r T (R+)- As for Q, it stems from the fact that 
divQ is quite smooth. 
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Therefore combining Corollary A.l, Lemma A.l and Lemma 2.2 allows to bound u(t) in 

• 2— - 

B Pt r r (W + ) in terms of the data and of the norm of VLT in L r (0, t; L p + ). In addition, because 
any function in L r (M + ; L p + ) may be approximated by smooth functions compactly supported 
in R + x R+, and because uq may be approximated by functions in W 2 ' p (R d h ) n Wq P (R+), 

■ 2-2 

the above formula guarantees that u is continuous in time, with values in B Pjr r (Wi_) (or in 

£p,r ? (R+) if divu = 0). 

In what follows, we assume that /j, = 1, which is not restrictive owing to the change of 
variables (2.13). Of course, when proving estimates, one may consider separately the three 
cases where only one element of the triplet (uo,f,g) is nonzero, a consequence of the fact 
that (2.1) is linear. 

Step 1. Case uq = and / = 0. Then the formula for u d given by Theorem 2.1 reduces to 



u d = Pg + U I e (t - T)A e a (Gk) dr, 
Jo 



/o 
and using the algebraic relations provided by Lemma 2.1 thus yields 

V 2 h u d = SUV h g + U I e^-^ A V 2 h e a (Gk) dr, 
Jo 

V h d d u d = (I- U)V h g +(I-U) [ e {t -^ A V h \D h \e a (Gk) dr, 

Jo 

dju d = (d d + (U-I)\D h \)g + f e^ A d d \D h \e a (Gk) dr + (I-U) [ e^ A A h e a (Gk) dr. 

Jo Jo 

The important fact is that all the terms corresponding initially to Pg may be written AVg 
where A stands for some 0-th order operator for which Lemma 2.2 applies. A similar ob- 
servation holds for the terms with the time integral so that applying Lemma A.l eventually 
yields 

(2-26) II W|| L r (L P) < C{\\Vg\\ L r T{LP+) + \\Gk\\ Lm) ). 

At this point, we use Lemma 2.2 to bound the right-hand side by ||(Vg, dtQ)\\ir <u>\, and we 
thus get 

(2-27) l|W|| L . (L , +) < C\\(Vg,d t Q)\\ L r T{LP+) . 

It is clear that V 2 u h also satisfies (2.27): indeed (2.10) gives 

(2.28) VV = r f e (t - T)A V 2 Se a {Gk) dr - SV 2 u d . 

Jo 

Let us now concentrate on the pressure. Keeping in mind (2.20) and (2.23), we may write 
d d U = (8 d - \D h \)g - S ■ Ud t Q h -(I- U)d t Q d - UGk 

-(r(\D h \-d d )\D h \ + UA) f e (t ~ T ^ A e a (Gk)dT. 

Jo 

Therefore, combining Lemmas 2.2 and A.l gives 

(2-29) \\dM L r T(LP+) < C\\(dtQ, Vg)\\ L r T{LP+ y 

Finally, because 

(2.30) V h U = -SGk - Sd d U, 
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it is clear that V/jII also satisfies (2.29). 

Step 2. Case / = and g = 0. With no loss of generality, one may assume that uq € 
Wl' p {R d + ) fl VF 2 ' P (IR+). From Theorem 2.1 and Lemma 2.1, we readily get 

V£tt d = t/e tA V2 ea (y d (no)), 
V fe a d u d = (/ - ?7)e* A V^|^|e a (y d (« )), 

9jn d = re* A ^|D^|e a (F d (u )) + r(I - U)e tA e a (A h V d (u )). 
Therefore, combining Corollary A.l and Lemma 2.2, 



IWH^) < c(|| ea (V 2 fc ^(« ))|| 2 



+ \\ea(V h \D h \V d (u ))\\._ hRd) + ||e.(a d |r> fc |^(ti ))||^_| CR(J) 

< C||uo|| .2-2 

Note that in order to bound the last term, we used the fact that because V d uo is null at the 
boundary, we have 

\D h \d d e a (V d (u )) = e s (d d \D h \V d (u )) € B~f(R d ). 
Owing to (2.28), \7 2 u h satisfies the same inequality. Finally, 

d d U = -(r(\D h \-d d )\D h \ + UA)e tA e a (V d u ) 
hence, according to Lemma 2.2 and Corollary A.l, 

ll^n|| L ^ (LP+) < C\\V 2 e tA e a (V d u )\\ L r T{LP) 

< C||l*o|| . 2 -2 

Bp,r T (IR+) 

Of course, (2.30) implies that V^LT has the same bound. 

Step 3. Case uo = and g = . As in the previous steps, owing to 

■ I e ( *- T)A V 2 e a (M/)dr - SVV, 
Jo 
V h U = SGf - Sd d U, 



V V = r 



it suffices to bound V 2 u d and d d U. The formulae for the second spatial derivatives of u a 
now read 

V 2 h u d = U f e^ A V 2 h e a (Nf)dT, 
Jo 

V h d d u d = (I-U) [ e^-^ A V h \D h \e a (Nf)dT, 
Jo 

dju d = r f eW^d d \D h \e a (Nf)dr + (I-U) f e^ A A h e a (Nf) dr. 
Jo Jo 

Therefore applying Lemmas 2.2 and A.l, 



IV 



2„.d 



v a hi(L> + )<C\\e a (Nf)\\ L r T{LP+) 

< C\\f\\ L r T{L v + y 
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For the pressure, we have 

d d U -f d = -UNf - (r(\D h \ - d d )\D h \ + UA) [ e^ A e a {N f) dr, 

Jo 

therefore, using once again Lemmas 2.2 and A.l, we obtain 

RII - f d \\ L r T (Ll) < C(\\UNf\\ L r T{LP+) + \\e a (Nf)\\ L r T{LP) ) 
< C\\f\\ L r T{LP+) . 

Step 4- Estimates for Vu. The starting point is the following classical Gagliardo-Nirenberg 
inequality on R : 

(2.31) \\Vz\\ Tmmd ,<C\\z\\ 1 - e _z d \W 2 z\ 

B p>r r (R ) 



| Lm(R d) ^<-IFIL 2 _2 II v z\\ LP(Rd) 



with 9 e (0, 11. m > p, < 1 - 2 + 20 < d d d. = d _ x + 2 _ 2£. 

If 9 € (0, 1) then this inequality may be easily proved by decomposing u into low and high 
frequencies by means of an homogeneous Littlewood-Paley decomposition (see e.g. [5] Chap. 
2 for the proof of similar inequalities). The case 9 = 1 corresponds to the classical Sobolev 
inequality. We omit the proof as it is standard. 

We claim that this inequality extends to the half-space setting if considering functions 
u S VFq' p (1R_|_) H W 2,P (M. + ). In effect, we observe that for such functions we have the following 
identities: 

V h {e a u) = e a (V h u) and d d (e a u) = e s (d d u). 

As V/j(e a n) also has null trace at (9R + , one can thus write 

V 2 h (e a u) = e a (y 2 h u) and d d V h (e a u) = e s (d d V h u). 

Even though d d {e a u) need not be zero at <9M + , it is symmetric with respect to the vertical 
variable, whence 

dj(e a u) = e a (dju). 

Applying (2.31) to z = e a u, and the above relations for the second order derivatives, we thus 
gather 



IV-uIIl" 1 < ||V(e u) 



\L m (R d ) 



< n\\p iiW 1 ^ \\xi 2 (p ii\\\ e 

- C||eatl|l ^-? C R-) 11 {eaun ^ d y 

< rWnii 1 - \\x7 2 i,\\ 9 

— ^ y ll"'ll 9 2 || * "'llrP* 

-Bp,r r (R+) 



Hence, taking the L q norm with respect to time of both sides (with q = r/9), we discover 
that 

Bounding the right-hand side according to the previous steps leads to the desired estimate 
of ||Vn|| L9(Lm) in (2.24). □ 

In order to solve System (1.4) for more general data, it will be suitable to extend the above 
estimates to the case where the index of regularity of no is not related to r. This motivates 
the following statement: 
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Proposition 2.2. Let 1 < p, r < oo and < s < 2. Let (uo,f,g) satisfy the compatibility 
conditions of Proposition 2.1, and be such that 



•'" >^ d \ t a (f XI n afll (= 77Y1U ,-TPA vrith rv = f 1 - £ - - 

2 r 



n € B s p , r (R d + ), t a (f, Vg, d t Q) € L r (R + ; L p + ) with a = 1 



Then System (2.1) has a unique solution (u, VII) with £ a (<9 t u, V 2 u, VII) € L r (R + ;Z,^) and 
for all T > 0, 



(2.32) |TOu, V 2 u,Vn)||^ (L , +) < C^||«o|lB| ir (R-) + ll* a (/> Vg,d t Q)\\ L r T{LP+) )- 

Furthermore, the following properties hold true: 

(1) For any couple (^2,^2) so that 

d d 2 2 

s < H <2-\ , P2>P, r 2 >r 

p pi T2 r 

we have t^Vu £ L r2 (R + ; L p + 2 ) with ft = \~ 4^ + 4f>~ h~ T^ and 

ll^v«|| 1?(L » ) <c(||iio|| i+ *. +\\nf,VgAQ)\\ L r T{Ll) } 

(2) For any couple (p3,r 3 ) so that 



d d 2 2 

s< <2H , P3>p, r 3 >r, 

p p3 r 3 r 



we have fu € L r *{R + ; L 1 ?) with 7 = -J- + |- - § - ± and 



2p3 2p 2 r3 ' 



II* 7 «IIl?(^3)<c(||mo|| | + ^ +||t a (/,V 5 ,^Q)|| L r (L p ) ^ 

Proof. Let us first assume that only g is nonzero. Then we start with the formula 

t a V 2 h u d = SUt a V h g + Ut a I e {t - r)A V 2 h e a (Gk) dr. 

Jo 

From Lemmas 2.2 and A. 2, we immediately infer that, if ar' < 1, 

\\t a V 2 h u d \\ L r (LP+) < C{\\t a V h g\\ L r (LP+) + \\t a e a (Gk)\\ LnLP{Rd)) ), 



whence 



(2-33) \\t a Vlu d \\ L r {LP+) < C\\t a (Vg,d t Q)\\ Ll{LP+) . 

Similarly, as 

t a V h d d u d = {I- U)t a V h g +(I- U)t a [ e^-^ A V h \D h \e a (Gk) dr, 

Jo 



and 



t a dju d = (d d + (U-I)\D h \)t a g + t a [ e^ A d d \D h \e a (Gk)dr 

Jo 

+{I-U)t a [ e^ A A h e a (Gk)dr, 
Jo 

it is clear that \\t a V 2 u d \\ L r / lP \ is bounded by the right-hand side of (2.33). Because 

t a V 2 u h = rt a [ e (t - T ^ A V 2 Se a (Gk) dr - t a SV 2 u d , 
Jo 
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the same inequality holds true for t a X7 2 u h . 

In order to bound the pressure, we use the fact that 

t a d d U = t a (d d - \D h \)g - t a S ■ Ud t Q h -(I- U)t a d t Q d - t a UGk 

-(r(\D h \-d d )\D h \ + UA)t a [ e^-^ A e a (Gk)dr. 

Jo 

Note that the terms in the right-hand side may be handled by means of Lemmas 2.2 and A. 2, 

exactly as we did for t a V\u . Therefore, we have 

\\edMu A L\) < c\\t a (v g ,d t Q)\\ L r {Ll) . 

Owing to (2.30), it is clear that t Q V/j!I satisfies the same inequality. 

Let us now consider the case / = 0, g = and u € W 1,P (M+) n W 2 ' P (IR+) (with no loss 
of generality). As usual, because one may go from u d to u h through 

u h = re tA e a (V h u ) - Su d , 

we concentrate on /r V 2 u . We start with the formula 

t a V 2 h u d (t) = Ut a e tA e a (V 2 h V d u ) = Ut a e tA V 2 h V d e a (u ), 

which, in view of Lemmas 2.2 and A. 5 ensures that 

s 1 



II^VMir^Pj < C\\V d e a (u )\\ss ARd) < C\\u \\^ pr{Rd+) with a = 1 - - - -• 

Similarly, we have 

t a V h d d u d = {I- U)t a e tA V h \D h \e a (V d u ), 
hence 

(2-34) \\t a V h d d u d \\ L r {LP+) < C\\u \\ BSpARd+y 

Finally, t a d\u d satisfies 

t*d 2 u d = rt a e tA e s (d d \D h \V d u ) + (I - U)t a e tA e a (A h V d u ) 
= rt a e tA d d \D h \e a (V d u ) + (/ - U)t a e tA e a (A h V d u ) 

because e s {d d V d uo) = d d e a (V d uo) owing to the fact that V d uo vanishes on 9M + . Hence 
t a d 2 d u d satisfies (2.34), too, and we conclude that 

(2.35) irvV|| W) <C||«o||^( K £)- 

Bounding VII is strictly analogous. 

In order to prove the estimate for t a V 2 u in the case g = and uq = 0, we use the fact 
that 

ev\u d = Ut a f e {t -^ A V 2 h e a (Nf) dr, t a V h d d u d = (I-U) t a f e^ A V h \D h \e a (Nf) dr 
Jo Jo 

and t a dju d = r I e {t - T ^ A d d \D h \e a (N f) dr + (I - U) t a f e^~^ A A h e a (N f) dr. 
Jo Jo 

Then combining Lemmas 2.2 and A. 2 readily gives 

rwii L , (LP+) < c\\t a f\\ L r {LP+) . 

Bounding t a V 2 Uh and t a VII works the same. 
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Let us finally go to the proof of estimates for t^Vn and fu. By virtue of (2.25) and of 
the definition of B (see the appendix), we have 

t?Vu{t) = r^Ve fA e fl (ji ) + t^Be a (f - VII) ' 

Therefore, applying Lemmas A. 3 and A. 5 yields 

\\iPVu\\l?(l'») < C[\\e a (u )\\^^ (Rd) + \\t a e a (f - VII)|| L[(LP(Rd)) J , 

whenever p 2 >p, r 2 > r, s 2 = ^ - f + s, 

d(\ 1 \ 1 1 1 1 a 2 1 

2\p p 2 J 2 r r 2 2 2 r 2 

and 

d d 2 2 

s < H <2-\ 

p p2 r 2 r 

Combining with the fact that e a is continuous on functions of Bp^ r . 2 (R + ) with null trace at 
the boundary, and with (2.32), we get 

\\^Al^{l^) < c (\\ u o\\B s P l r2 (Ri) + W ta ^^9,dtQ)\\ L r t{L v +) 
Finally, in order to bound t^u, we use the formula 

f u {t) = r^e tA e a (u ) + t^Ce a (f - VII) 
Applying Lemmas A. 4 and A. 5 yields 

W u \\L r t s (L p + s ) ^ C {\\ e a(uo)\\B s p l r3 (R d ) + W^^if ~ VII)|| l p (L r (K d)), 
with P3,>P, r 3 >r, s 3 = £ - ^ + s, 

d(\ 1\ 1 1 s 3 1 

2 \p P3/ r r 3 2 r 3 

and 

d d 2 2 

s < <2H 

P Pz r 3 r 

Combining with the fact that e a is continuous on functions of -Bi^ r3 (IR + ) with null trace, 
and with (2.32), we get 

This completes the proof of the proposition. □ 

3. Existence of smooth solutions 

As a first step for proving Theorem 1.1, we here establish the global existence of strong 
solutions for (1.4) in the case of a globally Lipschitz bounded density. As for the velocity, we 
assume that it has slightly sub-critical regularity. Here is our statement (recall that the space 
X p ' r has been defined just above (1.6)): 
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Theorem 3.1. Let a G W^°°(R%) and u G E^ t l + ' (R*) nB^'^) with p = 3^, 

r € (1, 00) and d < p < -35- ■ There exist two positive constants cq = co(r, d) arid ci = ci (r, d) 
so that if 

(3.1) ||^o||.-i + d < cq/j, and ||ao||i/>? < c\ 

B pr P (M + ) 

then (1.4) has a unique global solution (a, it, VII) with a G L~ (R+; W 1 > 0O (R^.)), 
(u, VII) G X p ' r n X p ' r and Vu G L 2r (R + ; L^' 1 ) n £ 9 (M+; L^°) with g 



d+(* -1)P 
In addition, there exist C\ = C\ (r, d) and C2 = C2 (r, d) so that 

(3.2) ||(u,VII)||jo>.r+M 1 "£||Vu|| ir < Cim 1- - ll^oll _i + d 

1 1 . 1_ 1_ . 1 1 

(3.3) \\(u,VU)\\ X p,r+H -^(\\Vu\\ L 2r (L1) +fi1 2r ||Vlt|| L «( L oo)) < C 2 U ~r\\u \\ . _ 1+ d , 

+ + Bp r P (K+) 

where a satisfies - = — + 2 -^- 

Proof. The general strategy is the same as in [10]: we set (o°,u°, VII ) = (0,0,0) and solve 
inductively the following linear system: 

( d t a n+l + u n ■ Va n+1 = 0, 

d t u n+1 - u.Au n+1 + Vn n+1 = F n 

(3.4) \ div u n+1 = 0, 

ii n+1 \ a — n 

, (a n+1 ,J +1 )\ t=0 = (ao,u ), 

with F n = f a n+1 (/iAn™ - VIP) - u n • Vn™. 

Both the global existence and uniqueness of a solution to (3.4) in the spaces given in 
Theorem 3.1 follows from basic results for the transport equation (given that Vu n is indeed 
in Lj oc (W + ; Ljf ) ) and from the solution formula for the evolutionary Stokes system in the half- 
space. So, as a first step, we focus on the proof of uniform estimates in L^ C (R + ; H /1,oc (IR + )) 
for the density, and in X^ r n XP< r for (u, VII). 

Step 1. Uniform estimates. It is obvious that 

(3.5) l|o n (*)lk~ = ||oo||l~ for all nGN and t G M + , 
and that 



(3.6) ||Va" +1 (t)|| L¥ < e / ° l|Vun|li + dT ||Voo|| 
We shall prove inductively that for all nGN, 

(3.7) ||K,Vn n )|| xP ,r + u 1 -^||Vn™|| dr <C^-r 

t t 2r I t 2r— 1 \ 



, ,Uo\\ . _ 1+ d 



(3.8) ||(n",Vn«)|| x!; , r +/, 1 -^(||V^|| Lr(L « ) + a^-||V^|| L?(L?) ; 



< CU 1 r|| Uo || _. 

B~ 



for some large enough constant C. 



'LftLl- 1 ) Lf(^+ ) 



GLOBAL SOLUTIONS TO INHOMOGENEOUS NAVIER-STOKES EQUATIONS 19 

Applying Proposition 2.1 we see that the left-hand sides of (3.7) and (3.8) are bounded 
(up to some harmless factor) by /U 1- ")!^ || _ 1+ d +\\F n \\ L r, L v\ and ix 1 ~~\\uq\\ _ 1+ d + 

\\F n \\ , p., respectively. Now, by virtue of Holder inequality (here we use that p = -^^ 
and the assumption p < -^- comes into play) , we have 

\F n \\ LULP+) < ||a n+1 || Lr(L?) ^||An"|| Lr(LP+) + ||Vn™|| Lr(LP+) ) + ||n"|| r2 * ||Vu* 

and 

lim^) < IK +1 || Lr(L oo)(>||Au"|| (LP - } + ||vn"|| ~ } ) + \\u n \\ * ||v«»||^. (L?) . 

Therefore, using the Sobolev embedding 

(3.9) W ( M +W L ~( M +)> 

we arrive at 

||(u n+1 ,vn n+1 )|| x p,.+^ 1 -27||vu n+1 || dr <cf/i 1_ 'll«oll -i+4 

I II 1+1 II / || A "-II I II V I I |7t|l \ i ||V7 "-l|2 \ 

+l|a llL t °°(L-)(Ml! An IIlj-(l^) + II vn \\ L r {L p +) ) + \\\/u || ^ J, 

L t r (L + ) 



H, 



" +1 I vn"+ 1 )|| x p + ^ 1 -i(||Vn" +1 || Lr(L « ) + ^-^||v^+ 1 || L?(L?) ) 

< C^-tWvoW 1+ * + \\a n+1 \\ LriL ^{n\\Au n \\ Lrt{L ^ + IIVIFH^) 

p,r *• -\- ' 

I llY7<,, r MI IIV7„i n ll 

+ VII dr Vtl \\ T 2r( T a\\. 

Now, the induction hypotheses (3.5), (3.7) and (3.8) thus imply 

||(u n+1 ,vn n+1 )|| x p,r + / u 1 -27||Vu n+1 || 

< C(l+ ||o ||i,°f + / t/~ 1 ||u || ._ 1+ d )^ 1_ -||u || __! 



B p>r (K_i_) -Bp.r 0&+) 



K+1)Vnn+1)|| ^ + ^__L (||Vn „ +1|| I^ ||Vuf , 



| x |?,r-t-M «-g|V"U llif(L«) "^ y " II v ^ + 1 1 L« (£,£>)) 



< C(l + llaolUs? +M 1 1 1 ""o 1 1 -i+d W HKH _ 1+ d 

B P , r P (K) B- r +P (K + ) 

Therefore, if c$ and c\ in (3.1) are small enough then we get (3.7) and (3.8) at rank n+ 1. 

Step 2. Convergence of the sequence. Let p be some real number in (^, 2 _i ) such that 
in addition 2 p < p < p. We are going to show that (a n ) n6 N and (u n , Vn n ) ng pj are Cauchy 
sequences in C&([0,T] xR + ) and Xj! r , respectively, for all T > 0. Of course, by interpolation, 
we easily find out that the bounds for (a n , VII ra ) ne N in X^ r are the same as in Xj! r n X^' r . 
As regards (a ra ) n eN, we use the fact that 

d t &i n + u n -V&i n = -5u n ~ 1 -Va n with 6a n = f a n+1 - a n and Su n d = u n+1 - u n . 



It would be natural to take p = p but we do not know how to handle the case r > 2 with this value of p. 
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Hence, using standard estimates for the transport equation, we get for all positive T, 

cT 



\Sa n (T)\\ L ~ < I \\Su n - 1 \\ L ^\\S/a n \\ L ^ dt. 

•J u 



Now, arguing exactly as in the proof of (2.31), we get the following Gagliardo-Nirenberg 
inequality 3 for functions z vanishing at d R , : 



(3.10) 
Hence 



\z\\l?° < C\\z\ 



«p,r r ( R +) 



|v «|| L p 



with 9 = I p ) — ■ 

2/ p 



(T)\\l?< / ||Va n || L oo||<5u 



n—l I 



2cn-l||l- 



\V z 5u 



(It 



/0 ^p./l^+J + 

Taking advantage of Young's inequality we thus get for all positive e, 

f-T 2p 

(3.11) \\5a n 



<r~ <rllV7 2 A, n— 1 II i fi / IIV7^, n— 111 ( 2 f>- d ) r \\X,,n— 1 II 

L°°(0,TxRl) - £ H V °" I ILL (L?. ) "+~ ° e / II Va HL°f II " IU2- 



St„ r (Ri) 



dt 



Next, we use the fact that 



d t 5u n - fiA6u n + V<JP = 5F T ' 
dw6u n = 



u n \ t= o = and u n 



tn 







with 



def 



5F n = L (fiAu n - VU n )6a n + a n {^8u n - 1 - Vdp- 1 ) - u n ■ V<5u n_1 - 5u n ~ l ■ V/' 1 
Applying Proposition 2.1, we see that for some constant Cq = Co(p,d), 

(3.12) \\(5u n ,Vm n )\\ x p,r+fi 1 -^\\V&i n \\ L 2r {Lm) <C \\5F n \\ Lr{L p ) with 



d _d _ t 1 
m p r 



Therefore, in order to prove the convergence of the sequence, it is only a matter of getting 
suitable estimates in L r (0, T]L P + ) for the terms in 5F n . We easily get 



\\SF n \\ L r (Lt) < ||& n || L cc (0l TxRi)ll^« B - VH^ 



+ \\a 



ILL(L^) 



in— 1 1 



^II^A&i™- 1 -Vol I 

-Hltt'll dr \\V5u n ~ l \\ T 2r (Tm ^ + \\5u n ~ i -\\ T 2r ITm *\\\Vu 
"L, T (L, + ) II UL, T (L, )\\ 



lL°°(0,TxR^.) 

*"|| dr I'™-"" 1 

with 4 m* such that Sr = I - 2 + ~ 



n—l I 



dr 



m' p 

At this point, one may combine (3.9), the following Sobolev embedding 



W, 



l,m/m>d 







L 



m* (mid 



and the conservation of the L°° norm of a n . We end up with 

W n \\ L r T{Ll) < (ll'fa n || L o O(0)TxR , ) ||(^, VH")|| xr + Haoll^lK&t"- 1 , Vdl- 1 ) 



,n-li 



+(||Vu n - 1 || 2 ^ +||Vn"|| ^ j Ll )\\V6u 

L r f(L + ) Lrf(L + ) 



n—l I 



\I#(L%))- 



It suffices to apply the corresponding inequality in E d to function e a (z) 
Here p < 2 t-i comes into play. 
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Keeping in mind the estimates that have been established in the previous step and assuming 
that the smallness condition over ao and uq in (3.1) is satisfied with small enough constants 
Co and ci, we eventually get 

C o\\SF n \\ L r T{ Ll) < \{\\{5a n -\vm n - l )\\ xp r +^\\V5a n - l \\ Ll r {LV ) 

+^ 1 ^||n || i .,_a ( ^ )ni . ! .2 (R:i) ||& n || L o O(0 , rxR:[) . 



p,r 



Now, plugging (3.11) with suitably small e in the above inequality, and resuming to (3.12) 
yields 



cT 



||(^,V<fI")|| xr +^ 1 -^||V ( 5u"|| Lr(L!F) < i (||(^- 1 ,V^ 

V^IIV^- 1 !!^™) +C I' A(t)\\(Su n -\ViR n - l )\\ x p,rdt 
+ Jo * 

where A is a continuous function of t that depends only on the norm of the data. Summing 

up over n > 1 and remembering that (u°, VIT ) = 0, we thus get 

£(||(& n ,V<IP)|| X £r +/x 1 -^||V<5u n || L 2 r(L!?) ) < ii^ 1 , vn 1 )!!^ + / i 1 -^||Vn 1 || Lir(L!r) 

n>l 

+4C f A(t)||(u 1 ,Vn 1 )|| X p. I .dt + 4C f A(t) S2\\(5u n ,Vm n )\\ x p,rdt. 

Jo « Jo „>! 

Hence, using Gronwall lemma implies that 

^2(\\(5u n ,vm n )\\ xpr + ^\\v5a n \\ Ll r {LV ) < (ll^.vn 1 )!!^ 

n>l ' ^ 

V-*l|Vu 1 ||^ (L?) +4Cjr T A(t)||(tt 1 > Vn 1 )|| Af , r dt)e 4Cf J?' A W*- 

This obviously entails that (u", VII n ) n gK is a Cauchy sequence in Xj! r for all T > 0. Then 
resuming to (3.11) implies that (a n ) n eN is a Cauchy sequence in C&([0, T] xl + ) for all T > 0. 

jStep 3. -End o/ the proof of the theorem. Granted with the convergence result of the previous 
step, and the uniform bounds of the first step, it is not difficult to pass to the limit in 
(3.4): we conclude that the triplet (a, u, VII) with a = lim n ^ +00 a n , u = lim„_> +00 u n , 
and VII = lim„^ +00 VIT 1 , satisfies (1.4) in the sense of Definition 1.1. In addition, as 
(n™, Vn n ) ng fj is bounded in the space X p,r n X p ' r which possesses the Fatou property, one 
may conclude that (u.VII) € X p < r D X p ' r and 5 that (3.2), (3.3) are fulfilled. Similarly, the 
uniform bounds for a n allow to conclude that a € L^ C (IR_|_; VF 1,00 (]R + )). This completes the 
proof of Theorem 3.1 □ 

4. Proving the existence part of the main theorem 

This section is dedicated to the proof of the existence part of Theorem 1.1. It is mostly 
based on a priori estimates for smooth solutions -the same as in the previous section, and 
on compactness arguments. 



Rigorously speaking we do not get the time continuity for u, but it may be recovered afterward from 
Proposition 2.1 by observing that u satisfies an evolutionary Stokes equation with source term in L r (R + ; IP, l~l 

IP,) and initial data in Sp~~ (RJ) n B-" 7 (R+) . 
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Step 1. Constructing a sequence of smooth solutions. This is only a matter of smoothing out 
the data (ao,Uo) so as to apply Theorem 3.1. We proceed as follows: 

• Let x ^ C£°(E d ) with x(x) = 1 for \x\ < 1. We extend oq to oq on R d by symmetry 
then use convolution of xix/ 71 )^ with a sequence of nonnegative radially symmetric 
mollifiers, then restrict to the half-space. We get a sequence (ao)neN in 1F 1,00 (E + ) 
with the same lower and upper bounds as ao, and satisfying ag — > ao a. e. on E + . 

• — 1+- 

• As no € B P: r P (E + ), one may take the zero extension uq = eo(uo) of uq over the 

whole space and then set 

\j\<n 

It is obvious that each term Uq is smooth and 7^0 = 0. Furthermore, («o)neN 

— 1+- • — 1+- 

converges to no in i? P)r p (E + ) n i?~ r p (E+). Of course, one may find p in (d, fry) 

.— 1+- 
so that each Uq belongs to B pr P (E+). 

Step 2. Uniform estimates. Let us solve system (1.4) with regularized initial data (oq,Uq) 
according to Theorem 3.1. We get a global solution (a n ,u n , VIP) in L~ C (E+; VK 1 ' 00 (E+)) x 
(I^nP' r ) satisfying 

II _ni| 1 1 Ti 1 1 <- || || 

II" Hl°°(R + xRJ) — H a 0llL°°(R^.) - ll a 0|lL°°(Rj)' 

and also 

.l-^llV7„.mi . / /^..1-i 



-1 



(4.1) ||(n n ,Vn n )|| X p,,- + n i -^||Vu n || *• < <V" \\ Uo \ { d 

In addition, by following the computations leading to (3.8), it is not difficult to see that the 
assumption that p > d is not needed if it is only a matter of getting a control on the norm 
of the solution in X' p,r . Therefore we also have 



(4.2) hk\vip)|| xp ~, <cy-i 



no|| d_ 1 

BE (Ri) 



p,r 



Of course, if p > d, then we also have a bound for V«" in L 9 (E + ; L5f). 

Step 3. The proof of convergence. Owing to the low regularity of ao, it is not clear that one 
may still use stability estimates in order to prove the convergence of the sequence defined 
in the previous step. In effect, as pointed out in the previous section, there is a loss of one 
derivative in the stability estimates for the density. Therefore, we shall use compactness 
arguments instead, borrowed from [14]. For completeness, we outline the proof here. 

According to the previous step, (dtu n ) n £^ is uniformly bounded in L r (E + ; L p + ). Combining 
with (4.1), (4.2), Ascoli-Arzela Theorem and compact embeddings in Besov spaces, we con- 
clude that there exists a subsequence, of (a n , u n , VII n ) ne N (still denoted by (a ra , n n , VII n ) ne N ) 
and some (a,u, VII) with a £ L°°(E+ x E|), 

dr ^ 

Vu £ L 2r (E+; L^- 1 ) and d t u, V 2 n, VII G L r (E + ; L\ n L\ ) 
(and also Vu G L q (R + ;L c £) if p > d), such that 

a n ^a weak * in L°°(E+ x«f), 

V 2 u n ^V 2 u and VIT 1 -»■ Vn weakly in L r (E+; L p + ), 



( 43 ) 
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with in addition for all small enough rj > 0, 

u n ^u strongly in L 2r oc (R + ;L^~ V (R d + )), 

' — n 

Vu n ^Vu strongly in L^(M+; L^ 1 '(R*)). 

By construction, (a n ,u n , VII n ) satisfies 

a n (d t (/) + u n -V4>)dxdt+ I <f){0,x)a%(x)dx = 0, 



(4.5) 



OO l> 

I <j) div u 11 dxdt = and 

JRi 



[ I \u n -d t <5> + ((u n ®u n ):V<5>) + (l + a n )(pVu n -Vn n )-<5>\dxdt 
Jo Jmf ^ > 

+ / v%-$(0,x)dx = 0, 

Jm.^ 

for all test functions 4>, $ given by Definition 1.1. 

Putting (4.3) and (4.4) together, it is easy to pass to the limit in all the terms of (4.5), 
except in a n (/j,Au n — VII n ). To handle that term, it suffices to show that a n — > a in 
L^ C (M + x IR+) for some m>r'. 

Now, it is easy to observe from the transport equation that 

d t (a n ) 2 + div(u n (a n ) 2 ) = 0, 

from which, (4.3) and (4.4), we deduce that 

(4.6) d t ~a? + div(ua2) = 0, 
where we denote by o 2 the weak * limit of ((a n ) 2 ) n ^. 

Thanks to (4.3), (4.4) and (4.5), there holds 

dta + div('ua) = 

in the sense of distributions. Moreover, as Vu € L (R + ; L^- 1 (R+)) and divu = 0, we infer 
by a mollifying argument as that in [17] that 

(4.7) d t a 2 + div(ua 2 ) = 0. 
Subtracting (4.7) from (4.6), we obtain 

(4.8) a t (^-a 2 )+div(n(^-a 2 ))=0, 

from which and Theorem II. 2 of [17] concerning the uniqueness of solutions to transport 
equations, we infer 

(a 2 - a 2 )(£, x) = a. e. iflj and f G R+ • 
Together with the fact that (a ra ) rae N is uniformly bounded in L°°(]R + x K + ), this implies that 

(4.9) a n -> a strongly in L^ C (K+ x R d ) for all m < oo. 

Granted with this new information, it is now easy to pass to the limit in (4.5). Therefore 
(a,n,VII) satisfies (1.7) and (1.8). Moreover, thanks to (4.1) and (4.2), there hold (1.10), 
(1.11) and (1.12). Besides, as (u, VII) satisfies (1.5) and the r.h.s. is in L r (lR + ; L^_nL^_), the 
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time continuity for u stems from Proposition 2.1. This completes the proof of the existence 
part of Theorem 1.1. 

5. More general data 

Until now, we assumed that p and r where interrelated through 

, , d 2 

5.1 -l + -=2-- 

p r 

It is natural to investigate whether the inhomogeneous Navier-Stokes equations may still be 

■ — 1+- 
solved with initial data (oo,«o) i n ^+ x &p,r P (R+) if (5-1) is not satisfied. 

The case where 1 < r < 2p/(3p — d) or, equivalently p < ^^ is not so interesting because, 

by embedding one may find some p* G (p, d) so that uq G <B„* ir P * (R ) and (5.1) is fulfilled 
by (p*,r). 

The case where r > 2p/(3p — d) is more involved and cannot be solved by taking advantage 
of embeddings. In order to explain how this may be overcome anyway, let us first focus on 
the toy case where u satisfies the basic heat equation 

d t u - Au = in M + xR (l 

— 1+— • —l+i 

with initial data u° G B p , r p (R )■ Then by using embedding in B~~ p (R ) for any p > p 

and r > r, we easily get (see Lemma A. 5 in the appendix) that 

(1) t a V 2 u G L r (R + ; LP(R d )) with a = f-^-±ifp> d/3, 

(2) t^Vu G L r2 (M + ; LP 2 (R d )) with /3 = 1 - ^ - ^ if p 2 > p, r 2 > r and p 2 > d/2, 

(3) ^w G L r 3(R+; L^(K d )) with 7 = I - ^j- - i if p 3 > p, r 3 > r and p 3 > d. 

As pointed out in Proposition 2.2, those properties are still true for the free solution to the 
Stokes system in the half-space. Keeping in mind that we want to apply those types of 
estimates to System (1.5), we see that we need to be able to handle also the Stokes system 
with some source term / satisfying t a f G L r (R + ; L p + ). Still in the simpler case of the heat 
equation: 

d t v -Av = f in R + x R d with t a f G L r (R + ; L p + ), 

it has been observed by the second author and collaborators in [14] (see also the Appendix) 
that if ar' < 1 then 

(1) t a V 2 veL r (R+;LP(R d )), 

(2) t?Vv G L^{R+;I^{R d )) with (3 = a+ f (± - £) - \ + \ - ± if p 2 > p, r 2 > r and 

d d_ < 1 1 _2 2 

p p2 ' r2 r ' 

(3) £^ e L^R+j L^(R d )) with 7 = a + f (± - -M - 1 + I - i- if p 3 > p, r 3 > r and 



2 V p p 3 / r r 3 



p 3 > d and jj - ^ < 2 + ^ 



In the case we are interested in, owing to the presence of u ■ Vu in (1.5) and to Holder 
inequality, the following supplementary relations have to be fulfilled: 

/ x 111111 

(5.2) - = 1 , - = 1 and a = (3 + 7. 

p p 2 p 3 r r 2 r 3 

Under the first two conditions, if se set 

,- „x def 3 d 1 def 1 d 1 def 1 d 1 

2 2p r 2p 2 ^2 2 2p 3 r 3 
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then the relationships above between a, ft and 7 are satisfied. Let us emphasize that if 
(p,r) with d/3 < p < d has been chosen so that r > 2p/(3p — d) (which is equivalent to 
< a < 1 — 1/r), then one may take any (^2,^2) such that P2 > p, P2 > d/2, T2 > r and 

, N d 2 2 d 2 

(5.4 --1 + < — <2 

p r r 2 P2 r 2 

The assumption on (p, r) ensures that such a couple exists. This motivates the following 
statement: 

Theorem 5.1. Assume that r G (l,oo) and p G (d/3, d) satisfy r > 2p/(3p — d). Then for 

any data Mo € B PjT v (M + ) and ao G £+ fulfilling the smallness condition (1.9), System (1.4) 
has a global solution (a,u) with \\a(i)\\L°° = ||ao||L°° f° r all t G R+, and 

(5.5) i a (3 t M,VVvn)eL r (l + ;L p + ), t^Vu G L r2 (R+; L^ 2 ), r« G L r « (K + ; L^ 3 ) 
with a,/?, 7 defined in (5.3), (p2>f2) satisfying (5.4), and (^3,^3) defined in (5.2). 

•2— — 

Furthermore, for any 1 < cr < 1 T^, , the fluctuation u — ul belongs to C(M. + ;B p>a ' y (R+)) 
where u/, stands for the "free solution" to the evolutionary Stokes system, namely 

(5.6) d t u L - u.Au L + VII L = 0, divuL = 0, itx,|t=o = u . 

Proof. For simplicity, we just treat the case \i = 1. The usual rescaling gives the result in 
the general case. 

We smooth out the initial velocity uq into a sequence («Q) nG N satisfying the assumptions 

.— 1+— 
of theorem 1.1: we take po so that — H = 2 — and require uo, n to be in B Pihr P0 (M + ) fl 

■- 1+ — j ■ -1+— rf 

i3 Pi r P0 (M + ), and to converge to u$ in 73 Pir p " (M + ). To construct wo.n; one niay first consider 

the zero extension uq = eo(uo) of % to 1 , then approximate it by compactly supported 

divergence free vector-fields by means of the stream function. 

Step 1. Uniform estimates. The corresponding solution (a n ,n n ,VII n ) satisfies in particular 

d t u n , V 2 u n , VIT G L r (R + ; V" n L\). 

Hence Holder inequality ensures that t a (d t u n ,V 2 u n ,VU n ) G L r (0,T;L p + ) for all T > 0. 
A similar argument ensures that t"Vu n and £ 7 u n are in L r2 (0,T;L^ 2 ) and L r3 (0,T; Zqf), 
respectively, for all T > 0. Now, because (a n ,u n , Vn ra ) satisfies 

a t u n - Au n + Vn n = F n = f a n (Au n - Vn n ) - u n • Vn n , 
Proposition 2.2 implies that 

Z n (t) d ^u n \\ L?{LP+3) + \\t^u n \\ L?{L?) + \\t a (d t u n , VV, VU n )\\ L r {LP+) 



< C {\\ u o\\ .-i+± . + ¥ * \\LULl))- 

Taking advantage of Holder inequality, of the relationship between (^2,^2) and (^3,^3), and 
of the conservation of ||a n (t)||i^>, we may write 

ll+ a Z? n ll <<■ ll„ II (\\+ a A„, n ll I ll+ a V7TT ?l ll \ 1 ll-/-7„, n ll \\+PY7„, n \\ 

¥ b \\L\iL\) < \\ a 0\\L™[\\t All || L r (L P) + ||t Vll \\L r t (L p + )) + ¥ u llLj 3 (L^)ll f Vu 1 1 1^ 2 (L* 2 ) • 

Therefore taking cq small enough in (1.9), we get that 



Z n (t)<C(\\u \\._ 1+si +2%(t)), 



B pr (1K_|_) 
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so that as long as ||uo|| ._ 1+ d < -^ji, we have 

B pr (R+) 

Z n (t)<2C\\u \\._ 1+ d d . 

B pr (K_|_) 

S'iep 2. Convergence. Hence sequence (a n ,u n , VII n ) ne ^ is bounded in the space (5.5). In 
order to complete the proof of existence, we have to establish convergence, up to extraction, 
to a solution (a, u, VII) of (1.4) in the desired functional space. For that, we first notice that 

Au n = r a (t a Au n ), 

Holder inequality guarantees that (Au n ) n( zfi is bounded in L° oc (R + ; If + ) for any a < r/(l + 
ar). Note that because a + 1/r < 1, one may take a > 1. Similarly, we have (dtu n ) n ^ 

and (Vn n ) ne pj bounded in Lf oc (R+; L p ,). Therefore, setting u n = u n — u\ where u\ stands 
for the solution to (5.6) with Uq instead of uq, we conclude that (u n ) n& ^ is bounded in 

C(R + ;Bl;J(R d + )). 

Similarly, writing that X7u n = t _/3 {t p Vu a ) and that u n = i~ 7 (t 7 M™), we get (W% eN 
and {u n ) n £fi bounded in L°" 2 (M+ ; L+ ) and L" 3 (K+ ; L^j 3 ) , respectively. As we may choose 
02 and 03 as close to (but smaller than) ^/(l + /3r2) and ^/(l + jr%) as we want, one may 
ensure that 

(5.7) — + — < 1. 

02 03 

Now, combining with the boundedness of (Au")„ e ^ in Lf (M.+ ; LPA and using Arzela-Ascoli 
theorem, we conclude that, up to extraction, sequence (a n ,u n , Vn n )„ g pj converges weakly to 
some triple (a,u,VU) with a G L°°(1R + xlj), 

f(a (J i,V\,Vn) G L r (R+;I^), t^Vn G L r ' 2 (R + ;LP 2 ), Fu G L r ' 3 (R+;L^). 
More precisely, we have 



a™ ->• a weak * in L°°(IR + x R%), 

V 2 u n ^V 2 u and Vn n ->■ Vn weakly in L a (R + ; L p + ), 
with in addition for all small enough n > 0, 

u n ^u strongly in Z,£(K+; ^""(O)' 
V^ -► Vu strongly in Lf 2 c (IR + ; Lf^K)). 

Because (5.7) is satisfied, passing to the limit in System (1.4) follows from the same arguments 
as in the proof of Theorem 1.1. That the constructed solution has all the properties listed in 
Theorem 5.1 is left to the reader. This completes the proof of existence. □ 

As in Theorem 1.1, assuming just critical regularity for the velocity does not seem to be 
enough to ensure uniqueness. At least, the constructed velocity u does not satisfy Vu G 
LJ oc (R + ; L+) so that we cannot resort to the Lagrangian approach. This motivates the 
following statement. 

Theorem 5.2. In addition to the hypotheses of Theorem 5.1, assume that uq belongs to 

■ — 1+- 
B~ r p (M_|_) for some hnite p > p. Then (1.4) has a global solution (a, u, Vn) fulhlling the 

properties of Theorem 5. 1 and, in addition, 

t a {d t uy 2 u,VW) G L r (R + ;Ll), t^Vu G L r2 (R+ ; L^ ) , f f u G L r ' A (R+ ; L^ ) 
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with a, f3 and 7 defined as previously in (5.3), 

11 - 11 r/ 11 - 11 

P2 P2 P P PZ P3 P P 

whenever p 2 and r 2 may be chosen so that (5.4) is fulfilled and 

, s 2d d d 2 2 

5.9 <1 + 

p p p 2 r 2 r 

If in addition p > d then there exists some positive 5 < a and s > r so that 

(5.10) t j VuGL s (R + ;L"), 

and the constructed solution is unique in its class of regularity. 

Remark 5.1. Despite the appearances, it is always possible to take p > d in the above 
statement. At first sight it seems not obvious because a necessary condition for having p > d 
in (5.9) is that p > 2r -i - However, by embedding, one may always find some p\ £ ( 2r ^ 1 , d) 

so that n G B pur Pl (R+) nB~ r Pl (R^) and thus replace p by p x . 

Proof. The general scheme for proving existence is exactly the same as in the previous state- 
ment. Therefore we remain at the level of a priori estimates. Let Z and Z be defined on R + 
by 

z(t) A M ||ru|| L?(L p 3) + ||^Vu|| L?(LP+2) + \\t a (d t u, v 2 u,vn)\\ LULP+) 

Z(t) t f \\Pu\\ L?{L?+3) + ||^Vn|| L?( ^ 2) + \\r{d t u,V 2 u,VYi)\\ L r {L v +y 

According to the proof of Theorem 5.1 and under Condition (1.9), we have for some constant 
C = C(p,r,d), 

(5.11) Z(t) <C||«o|| ,_ 1+ a d ■ 

B pr (K_j_) 

Next, keeping in mind our assumptions on the Lebesgue exponents p, p 2 , r 2 , p and p 2 , 
Proposition 2.2 ensures that for some constant C = C(p,p,r,d), 



Z{t)<C(\\u \\ _ 1+ d +\\t a a{Au-VU)\\ Trlr p, + \\t a u-Vu\\ rrlr ^). 

B- p (R d + ) " • 



~p,r 



Because ||<x(t)||i°° is constant during the evolution, under Condition (1.9), the first term of 
the right-hand side may be absorbed by the left-hand side. As for the last term, we use 
Holder inequality and the fact that 



111 j ! 1 1 

a = p+7, - = 1 and - = — H 

r r 2 r 3 p p 2 p 3 



We thus end up with 



Z(t) < c(\\uo\\._ 1+ $ + II^VnU^^^Tnll^a^)), 

p,r *• ~\~ ' 

whence, if c is small enough in (1.9), 

(5.12) Z(t) <C||«o|| _ 1+ d 

B- r ?(Ri) 

In order to prove (5.10), we first have to check whether one may take p > d, knowing that 
Conditions (5.4), (5.8) and (5.9) have to be fulfilled. This is in fact equivalent to p > 2r -i - 
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Assuming from now on that this condition is fulfilled, and taking p > d, we may use the 
following Gagliardo-Nirenberg inequality for all functions 6 of Lf f] W 1,P (M+) n IU 2 'P(M+): 

(5.13) l|Vn|| L oo<C||V 2 nf r? ||^ir^ with 6 = f l + F " ~/ 2 • 
Then, using Holder inequality and (5.12), we readily get 

l|t 5 vu|| L?(L?) < c\\evM{r {Ll) \Wu\\^ {L , Z) < ciitioii 1+ . < oo 

with 

(5.14) S d = 0a + (1 - (9) 7 and - = - + l ~ ° 



'.) 



s r T3 

Because r^ > r and 7 < a, it is obvious that s > r and <5 < a. This completes the proof of 
the second part of the statement. Proving uniqueness is postponed to the next section. □ 

6. Uniqueness 

This section is devoted to proving the uniqueness parts of Theorems 1.1 and 5.2. As in 
[11, 14], it strongly relies on the fact that for smooth enough solutions, one may use the 
Lagrangian formulation of (1.4), which turns out to be equivalent to (1.4). 

6.1. Lagrangian coordinates. Before going into the detailed proof of uniqueness, we here 
recall some basic facts concerning Lagrangian coordinates. Throughout, we are given some 
smooth enough solution (a,u, VII) to (1.4) (typically we assume that (u, VII) € X^ r n Aj,' r 
with (p,r) and p > d as in the statement of Theorem 1.1). Then we set 

b(t,y) d =a(t,X(t,y)), v(t,y) d ^ u(t,X(t,y)) and P(t,y) d = U(t, X(t,y)) 

where, for any y G M. , X(-,y) stands for the solution to the following ordinary differential 
equation on [0, T] : 

(6-1) dX{ ^ y) = u(t, X(t, y)), X(t, y)\ t= o = y. 

Therefore we have the following relation between the Eulerian coordinates x and the La- 
grangian coordinates y: 

(6.2) x = X(t,y)=y+ f v{ T ,y)d T . 

Jo 

Let Y(t, •) be the inverse mapping of X(t, •), then D x Y(t,x) = (D y X(t,y))~ l with x = 
X(t,y). Furthermore, if 

(6.3) / ||Vu(t)||i,«>dt<l 

Jo 

then one may write 

(6.4) D X Y = (Id +{D y X -Id))" 1 =^(-l) fe (/ D y v(r,y)dr 

k=o ^° 



As know usual, this inequality may be deduced from the classical one on R , taking advantage of the 
antisymmetric extension operator. 
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Setting A(t,y) = (D y X{t,y))~ l = D x Y(t, x) for x = X(t,y), one may prove (see the 
Appendix of [11]) that 

(6.5) V x u(t,x) = T A(t,y)V y v(t,y) and div x u(t,x) = div y (A(t,y)v(t,y)). 
By the chain rule, we also have 

(6.6) dw y (Av) = T A : V y v = D y v : A = Tr (D y v ■ A). 
As in [12], we denote 

V u = A ■ Vy, div u = div(-A-) and A u = div u V u . 

Note that for any t > 0, the solution of (1.4) obtained in Theorem 1.1 satisfies the smoothness 
assumption of Proposition 2 in [12], so that (b,v,VP) satisfies 



(6.7) 



( bt = in R+ xKf, 

R+ xK{, 
R+ x Rf, 



d t v-(l + b)(A u v-V u P) = in 
div u v = in 

v = on R+X3R+, 



-oil 



, (6,u)|t=o = (oo,«o) in 

which is the Lagrangian formulation of (1.4). 

6.2. Proving uniqueness : the "smooth case". Here we prove the uniqueness part of 
Theorem 1.1. Thanks to the usual rescaling (2.13), one may assume with no loss of generality 
that n = 1. Let (aj,Uj,Ilj),i = 1,2, be two solutions of (1.4) satisfying (1.10), (1.11) and 

(1-12). 

For i = 1,2, let X{ be the flow of ut (defined in (6.2)) and denote by (vi,Pi) the corre- 
sponding velocity and pressure in Lagrangian coordinates. Let 

Su= i v 2 -v 1 , SP = P 2 -Pi. 



Observing that b = ao, we see that (Su, V5P) satisfies 

(6.8) 



d t 6v -A6v + V5P = a {A6u - V5P) + 5f\ + 5f 2 = f SF in 



where 



(6.9) 



div dv = Sg = div 5R in 



Sfi d = (1 + a )[(Id - T A 2 )V5P - 5AVP X } with SA d = A 2 - A u 
Sf 2 d = (1 + a ) dW[(A 2 T A 2 - Id)V<5u + (A 2 T A 2 - A 1 T A 1 )Vv 1 ], 



Sg ^ (Id - T A 2 ) : V* - T 5A : Vv u 
SR d = (Id- A 2 )5v - SA Vl . 



Here and in what follows, we always denote by A the transpose matrix of A, and (Vti) . . = (ftu J ) . .. , 
and Du = T Vu = (djU 1 )^. .^ J . 
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As ~j5R d = and &>|t=o = 0, applying Proposition 2.1 with p = 3 _ 2 , and (q,m) = (2r, 2 _i ) 

2r-2/ 



or (g,m) = (r, 2^2) implies that 



-1- |((//or. \"- 



|V<5u|| Tr dr T +IIV&H ir dr Tr + ||(d t <5y,V 2 <5y,Vff , )|| I ,r (I P ) 



< C(J|a (A& - V5P)\\ L r {LP+) + \\(Sf lt Sf 2 , Vfc,$ffi)|| L r ( ^ 

Of course, the first term in the right-hand side may be absorbed by the left-hand side if the 
constant c\ is small enough in (1.9). So let us now bound the other terms. In what follows, 
we will use repeatedly that (see e.g [11]) 

ft ft 

(6.10) SA{t)=U DdvdT^(j2 Yl Ci'C*" 1 " 5 ') with Ci{t) d = Dvidr. 

Bounds for 6g . The definition of 6g implies that 

l|V&|| L r (L P) < ||V^ 2 ®V*|| L , (L P ) + ||(Id-^ 2 )^V 2 *|| L , (L P ) 

+ ||V&4 (g> Vvi||£r( L P) + II &4 <g) V 2 Wl|| L r( L P-). 

Therefore using Holder inequality, (6.3) and (6.10), we get 

l|V4r|| L r- (L P) < cn|V^ 2 || Lr ( L d)||V*|| ^ +ll(Id-^2)||i,~(L«')||V 2 &;|| L r (L p ) 

L r (L, ) 



+ ||Vd4|| L «,( L P ) ||Vul|| L r( L oo) + ||(SA|| ^ l|VV||Lr (L d) 



72 

dr II V 
F=Z J 1 

Remark that for 9 being determined by \ = - + ^~-, 



Lf{Lf—') ' 



d p 



I v ^2 II l°°(l*) < C / ||V 2 v 2 ||i«t ^' < C't 1 "^||V 2 t; 2 ||^ (LP) ||V 2 t; 2 



2„. nl-fl 

Lr(j*)' 



.1-^ 



||(Id-^ 2 )||L»(L») ^Ct 1 «||V«2|| L?(L oc), 

||VM|| Lr(LP+) < c^-^iv 2 *!!^^), 

||&4|| dr < Ct 1- ^ || V&;|| *■ , 

l|vV|| Lr(Lt) <c||v 2 ,j^ (LP+) ||vV||^ LP;) . 

Hence we obtain 

HV^IIi^iP) <7/i(i)(||V<HI r ^ +||V 2 ^|| l , (l p ) ), 

with 771 (£) — )■ as t goes to 0. 
Bounds for dfSR. First, we see that 

||^[(Id-A 2 )&]|| L r (L P ) <C(||V^®&|| L r (L P ) + ||(M-A 2 )^&|| L r (L P ) ). 

It is easy to check that 

||VU2 <8> Su\\l v (L p ) — C||VU2|| dr \\Sv\\ dr 

<C||VU 2 || dr \\Vdu\\ dr , 
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and 



\\(ld-A 2 )dtdu\\ L r {L p ) <C\\U-A 2 \\ Lr ^)\\dtdv\\ L r^ L p ) 



,1- 



<Ct ^\\Vv2\\ L <3 {L ^ ) \\d t du\\ L r {L P +) , 
which gives rise to 
\\dt[(U-A 2 )dv}\\ L r (L P ) <r]2(t)(\\V5u\\ *_ + \\d t Su\\ L u L P)) with limrj 2 (t) -- 

11 +' L 2r (L ) *^"0 

On the other hand, thanks to (6.10), we have 

||3t[&4vi]||i,r( L P) < Ci \\v\®V5v\\^( I p + } + \\&Ad t vi\\ L r( L Pj+ / |V&>| dt'\Vvi^\\vi\ 

where v\ y2 designates components of v\ or v 2 . 
Applying Holder and Sobolev inequalities gives 

\V\ % VSu\\ L r( L P) <C||Vl|| dr ||V&|| dr 

< CIIVVlll dr ||V<fc|| dr . 

Following the computations for bounding 5g, we also have 



L r t (L P +) 



\\SAdtVlWL2rtL*) <C\\5A\\ dr H^VlHir/id) 



<Ct 1 ~r\ 



rpAldtvA 1 ' 



\V5v\ 



ytVl ^ULl)\^^ L r {L P +) U "^^ 



and, for a so that — \- '^-j— = -, one has 

' a dr » ' 



/' 

JO 



\V5v\ dt'\Vvi t2 \\vi 



L r t(L p +) 



< C||V&|| dr IIVVl 2 <8> Vl\\ T r(Td\ 

< Ct^\\Vdv\\ dr (||Vfl2|| dr \\V\\\ dr 

_ y- e . 



X 



(||Vvi,2||l,2r( L c<)||vi|| _Ar_ 



As a consequence, we obtain 

II <9 t [M vi] II rwrp-, < r/ 3 (t)(||V(i;j| dr + ||V&>|| dr ) with lim77 3 (t) = 0. 

Bounds for 5f\ . We notice that 

5h = (1 + oo) ((Id - T A 2 )V5P - T MVPi) . 
Hence, thanks to (1.10), (1-11) and (1.12), we have 

WiWlulD < C(\\ T (Id-A 2 )VSP\\ L r {L P +) + \\ T 8AVP 1 \\ L r, L P ) ) 



< C(||Vt>2||, 1(T °°\\\V 5P\\ r r ( r P \ + HV&H dr || V P\ II r r I r d \) 



<C7^||V W2 || L?(L?) ||Vff'|| L r (LP+) +^|tVP 1 || LnLP+nL , ) ||V ( 5 U | 



m.K ) 



so that 



h\\ L i(Ll) < V4(t)(\\V5P\\ L r {LP+) + 



IV&I 



V t (Lf^)> 



with linnet) = 0. 
t— >o 
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Bounds for 5f 2 ■ We may write 

\W2WmLD < c{\\dw((A 2 T A 2 -Id)V5u)\\ L r {L P +) + || div((,4 2 T ,4 2 - ^i T -Ai)Vvi)|| L r (L P) 
Therefore Holder inequality and (6.10) imply that 

T 



Sf2\\ L r {L p +) <C(\\V(A 2 1 A 2 )\\ L oo {Ld+) \\VSu\\ ^ +||(^ 2 M 2 -Id)|| L oo (L oc ) ||V 2 (5 y || LKL P ) 



^Vv 1 \\^ iL ~ ) \\V(A 2 I A 2 -A 1 I A 1 )\\ L ~ (L P ) + \\V z v 1 \\ L r (L d ) \\A 2 1 A 2 -A 1 1 A 1 \\ 



dr 
ir-i - 



which leads to 

¥h\\Li{L\) <c(^(l|V 2 vi|| Lr(L d ) + ||V 2 v 2 || Lr(L d ) )+^(||Vvi|| L?(L oc ) 



+ I|Vv2||l«(l»)))(I|V&|| t ^ + ||V 2 &|| L r (L P )) 



<V5(t){\\V5u\\ or + ||V 2 (S;|| L r (L P )) with limr) 5 (t) = 0. 

Therefore, one may conclude that 

IIV&II *■ +IIV&II * +\\(d t 5v,V 2 5v,V5P)\\ L r (LP) 

where rj(t) = Xli=i Vi(t)i which goes to zero as t — )■ 0. This yields uniqueness on a small time 
interval. Then a standard continuation argument yields uniqueness on the whole interval 
[0, T]. The proof of Theorem 1.1 is now complete. 



6.3. Proving uniqueness : the "rough case". We now assume that we are given two 
solutions (ai,Ui,Tli),i = 1,2, satisfying the properties of the last part of Theorem 5.1. Of 
course, the difference (5u, VdP) between the two solutions in Lagrangian coordinates still 
satisfies System (6.8). In order to prove uniqueness, we shall derive suitable bounds for the 
following quantity: 

SZ(T) t f \\t a (d t du, V 2 dv, V5P)\\ L r T(LP+) + \\^VSv\\ L?{L r +2) + \\Pdu\\ L?{L i*y 

To start with , let us apply Proposition 2.2 with regularity exponent — 1 -\ — . We get for all 
positive T: 

(6.11) ^(T)<C(||a t Q (A&-V(5P)|| L r (L P ) + \\t a (5f\, 5f 2 ,V6g,d t 5R)\\ L r T{LP+) . 

Again, the first term of the right-hand side may be absorbed if the constant c\ is small 
enough in (1.9). 

Bounds for V6g . From the definition of 5g, we readily have 

irV^H^P) < \\t a VA 2 <g> V*|| L r (L P) + \\t a (U-A 2 ) V 2 5v\\ L r T(LP+) 

+ \\t a V5A ® Vui|| L r (l p ) + \\t a SA <g> vV||l5, (l p )• 

Using Holder inequality and (6.3), we get 

f T 

\\t a VA 2 ®V5v\\ L r (LP) <\\t a V5v\\ rruP *, / \\V 2 v 2 \\ T d dt, 

11 * \\l, t (l, + ) — 11 "L r T (L' J + ) J Q 11 ii^ + 
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def 
where p* = dp/(d — p) stands for the Lebesgue exponent in the critical Sobolev embedding 

(6.12) W llP (R^)^LP*. 
Therefore, remembering that a + 1/r < 1 and that L^_ ^-t> L p + n L p + , 

rvA 2 ®v&||^ ( ^ ) <crW-«||t«v 2 &||^ ( ^ ) rv 2 ^||^ ( ^ r ^- ) . 

Next, using the fact that 

(6.13) ||(Id-A 2 )(r)|| L oo <C / \\Vv2\\L~dt<CT 1 - k s- 5 \\t s Vv 2 \\ L s T{ L~), 
where s and 5 are defined in (5.14), we easily get 



|t Q (id-^ 2 ) ® v 2 *|| L ^ (LP+) < cr 1 -7- 5 ||^v V2 ||^(L-)||t Q v 2 



LJ(L5.)- 



I I Llrp y±J I 

In order to bound the third term of V&/, we notice from (6.10) that 
(6-14) ||VM|| Lm) < Ct l --r' a \\t a V 2 5v\\ L r T{LP+y 

Therefore, setting ^- = - — -, 

' ° m r s * 

< CT 1 -7-«||t«V 2 &;|| L r (L P ) ||t*V Vl || L| , (L =o ) . 
Finally, we have from (6.10) and Sobolev embedding (6.12), 

(6.i5) II^Hl^lO ^ c rl """ a IK Q v 2 *ll^(^ +) . 

Hence we obtain 

\\t a 5A® W|| L . (LP+) < cr 1 -?- a ||t a v 2 <j v || L5 , (L P ) ||« a v 2 t;ilLr ( ^ ) . 

Putting together the four above estimates, we conclude that 

+QV72£,,|| „ ( rT<\-\-a.\\.a.cr-,2„, ^-,2 



(6.16) \\t a V&j\\ L r T{LP+) < C||t Q V^||^ (L!h) ^-"«||r(W, V^JH^^p^pj 

+ T 1 -i- a ||t 5 (Vt; 1 ,Vt; 2 )|| w?) ). 

Bounds for dt5R. Recall that 

$ffi = -($A 2 )& + (Id -A 2 )d t 6u - (d t SA) v x - SAdm. 
Now, using the expression of dfA 2 and Holder inequality, we get 

\\t a (d t A 2 )5v\\ L r T(LP+) < C\\t^V 2 \\ L?{Ll2) \\t^Su\\ L r T S {L r +3 y 

Next, we have, using (6.13) 

\\t a (Id-A 2 )d t dv\\ L r T{LP+) < C\\Id-A 2 \\ L «>( L ~)\\t a dt&)\\ L r i ( I r + ) 

< cr 1 -i-*p*v^|| L5 ,( L =o)||t a ^&;|| L5 , (L P ) . 

In order to bound the third term of dfSR, we differentiate (6.10) with respect to time and 
easily find that 



\\t a (d t 5A)vi\\ L r {L P ) <c(\\t a vi®V5v\\ L r {L v ) + t a [ |V<5u|^'|Vvi, 2 ||vi| 
where v\ 2 designates components of v\ or v 2 . 



L* T {L* + ) 
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On the one hand, applying Holder inequality gives 

||£ Q Vl<g>V&||i-(Lp < C 'll t7v lllL^(L^ ) ||t /3 V(fc|| L?(L P2 ) . 

On the other hand, we have, by virtue of (6.12) 

t a / \V5u\dt'\Vvi 2 \\vi\ <C\\t a Vv 1 2®vi\\ Tr(T d ) \\V5u\\ rP *dt 

<CT 1 -7-«||t«V 2 &;b 5 , (L P ) ||t^ 1 || L?(L P3 ) ||^Vt;i 1 2|L ?(L a ) 

with i = i — J-. Given that p < d <p, we have P2 < a < p 2 , hence L™ ^-> L+ 2 n Zq 2 . 
Finally, arguing as for bounding the last term of V6g yields 

Putting all those estimates together, we conclude that 

(6.17) II W| LW) < r](T)5Z(T) 
for some function 77 going to at 0. 

Bounds for 5f\ . We notice that 

Sh = (1 + ao) ((Id - T A 2 )VSP - T MVPi) . 

Hence it suffices to follow the computations leading to the bounds for the second and fourth 
terms of V<5?: we just have to change V 2 <5u and V 2 t>i into V5P and VP\. We end up with 

(6.18) \\t a 5h\\ L r T(LP+) < c{T 1 -^- 5 \\t 6 Vv 2 \\ L s T(L ^ ) \\t a V5P\\ L r T{LP+) 

l T"! a ll+ Q V7D II ll+ a V72l,,|| | 

+ T r \\ t VPl\\ L r T(LP+r}L V +) \\t V 50\\ L r T{L P +) y 

Bounds for 5f 2 ■ As we may write 

II^IIl' (l" ) < c(\\t a V(A 2 T A 2 ) V5v\\ L r (l p } + \\t a (A 2 T A 2 - Id) V 2 &|| L 



+ ||t a V(^ 2 T ^2 - ^l^l) <8> Vvi|| L r (L P) + ||(^2 T ^2 - ^l T ^l)VV|| L ^ (L P)), 

one may repeat the computations leading to (6.16): this only a matter of replacing everywhere 
A\ and A 2 by A\A\ and A 2 T A 2 , respectively. We conclude that ||£ Q <5/2||i7 il v ) i s bounded 
by the right-hand side of (6.16). 

Conclusion. Plugging (6.16), (6.17), (6.18) and the above inequality in (6.11), we conclude 
that whenever both v\ and v 2 satisfy (6.3), we have 

5Z(T) < r){T)5Z(T) 

for some function r\ going to at 0. This implies uniqueness on a small time interval. Then 
a standard continuation argument yields uniqueness on the whole interval where the two 
solutions are given. 
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7. Remarks on the bounded domain case 

This section aims at extending partially our main theorem to the initial boundary value 
problem (1.4) in a C 2 bounded domain $7 of R . Before we present the main result, let us 
introduce a few notation. We set 

X q (n) d = { u £ (L q (n)) d \ divu = and u-n = on dfl }, 

where n stands for the unit normal exterior vector at dfl. Denoting by P q the projection 
operator from L q (£l) onto X q (Q), the Stokes operator on L q (£l) is the unbounded operator 
(see e. g. [13]) 

A q d = -P q A with domain D(A q ) d = W 2 > q {n) n W 1,9 (fi) n X q (Q). 
Definition 7.1. Let 1 < q < oo. For a £ (0, 1) and s G (1, oo), we se£ 



IU.II -^IU,II j- f f°° \\A-<xa c -tA„„,\\s dt\ s 

\\ u \\dY — \\ u \\li + I / ||* A?e 9 «||l9— I , 

where e~ q stands for the analytic semigroup generated by A q . We then define the inhomo- 
geneous fractional domain D^ s as the completion of D(A q ) under ||tt||rj a > s - 

Remark 7.1. Let a € (0,1) and 1 < q,s < oo. Let B q)S be the completion of C^°(Q) in 
B^ s (R d ). Then we have 

flg(n) n x q (n) ^ D a / q ^ sg(n) n x q (n). 

If moreover 2a < 1/q then the three sets coincide. 
The main result of this section reads: 



Theorem 7.1. Let Q be a bounded domain with C boundary. Let r € (1, oo), oq € L°°(0) 

l—-,r .] 

and uq € D A r ' with p = ^ r _ 2 • There exists a positive constant cq so that if 
(7.1) MlKllz,°°(fi) + 1 1 t*o 1 1 i-i, r < c 0/ u, 

*-* A 

then (1.4) has a global weak solution (a, u, VII) in the sense of Definition 1.1, which satisfies 

(7.2) 



^ 'Nl roor „ n i-i,, +\\(d t u,fiV 2 u,VU)\\ L r (M+ . LP{n)) 



+ it2r||Vli|| dr <Cu} r ||Uo|| 1 1 _, 

ll L2'-(M + ;L5F=T(n)) ~ ^ U D A "' 

for some sufficiently large positive constant C . 

Remark 7.2. Uniqueness would require our using Lagrangian coordinates, hence investigat- 
ing the evolutionary Stokes system in a bounded domain with non homogeneous divergence. 
We Jeave this interesting issue to a future work. 

The proof of the theorem mainly relies on the following result (see Theorem 3.2 of [9]): 

Proposition 7.1. Let Q, be a C 2 bounded domain of R and 1 < q,s < oo. Assume that 

u € D l ~ 1,s and f € L S (R+; L«(fi)). Then the system 



d t u - fiAu + VII = 


= / 


in 


div u = 




in 


u\an = 




on 



cfi, 
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with initial data uq has a unique solution (u, II) satisfying j Q II dx = and for all t > 0, 
A* 1-7 IM*)II i-i,. + \\(dtu, fiV 2 u,VU)\\ LtiLqin)) < C(/z 1- s||u || .i,, + ||/||L|(L«(n)))- 

D A q D A q 

Now we are in a position to prove Theorem 7.1. 

Proof of Theorem 7.1. We first solve System (1.4) with regularized data, according to e.g. 
[9]. We get a sequence (a n ,u n , VII n ) ne N of smooth solutions to (1.4). In particular, as 

d t u n - u.Au n + VIP = F n d = a n (jiAu n - VET) - u n ■ W\ 
Proposition 7.1 implies that 

Z n (t) d ^^\\u n (t)\\ ,_!,„ + ||ft« n ||L ? (LP) +Ml|Au n || LKL p) + ||Vir|| L r (LP) 



<C(^\\U \\ hr + \\F n \\ L r {Lp) ). 



D A 



By interpolation and embedding we have for all z € W 1,p (fl) the following Gagliardo- 

Nirenberg inequality: 

i. i 

IMI j < n\\?\\ 2 llrll 2 

Applying this inequality to z = V« n and taking advantage of Poincare inequality (here we 
use that u n vanishes at the boundary), we discover that 

i i 

||Vu n || *. <C\\u n \\\ 2 ||VV||», ov 

Therefore, using Remark 7.1, 

(7-3) IIV« B ll^ (Lafe(n)) ^ ^(^-Iholl^i, + \\F-\\ Ll(LPm ). 

Now, taking advantage of Holder inequality and of the Sobolev embedding 

1 ^ r rl 

w 3r - 2 (n) ^ L^(n), 

we may write 

II Z? n ll << ll^ n ll f,,IIA„, n ll l llT7TT n ll "\ 

ll-r \\Ll(LP(n)) S \\a ||£,°°(L«>(n)HM||AM ||l[(Lp(Q)) + II Vll ||L[(LP(n))J 

+ ||n n || dr ||Vu n || dr 

< !|aollL-(n)(/^l|Au n || Lr(L p (n)) + ||Vn n || LI(L p (n)) ) +C||Vn"||^ r ^^_ ( ^. 
Therefore taking Co small enough in (7.1), we get by using (7.3) that 

Zn(t) < CV-'IIUoLi-i,,. +li-<- 2 -r)Z%(t)), 

^ A 



so that as long as ||uo|| i_i r < i^nA 4 ) we have 



D 



j> 



Z n (t) < 2C/i 1 HKoll l-I r - 

Granted with this estimate, we can follow the lines of the proof of Theorem 1.1 to complete 
the proof of Theorem 7.1. □ 
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Appendix A 

Here we establish several L p — L q or maximal regularity type estimates involving the heat 
semigroup in the whole space, or the Stokes semigroup in a bounded C 2 domain fl. Although 
those estimates belong to the mathematical folklore (as a matter of fact the heat semigroup 
case in R has been treated in [14]), we did not find any reference where they are proved 
with this degree of generality. 

As in [14], the key to the proof of maximal regularity estimates is the following proposition 
(see e. g. Th. 2.34 of [5]) enabling us to characterize Besov spaces with negative indices by 
means of the heat semigroup. 

Proposition A.l. Let s be a negative real number and (p,r) G [l,oo] 2 . A constant C exists 
such that for all [i > 0, we have 



C ^ 2 ll/II.Bs r (Rd) < \\\\t 2 e^ f\\LP(R d )\\L^(R + ;f ) - C^ 2 llJ \lB^ r (R d )- 



We shall often use the following consequence of the above proposition: 

Corollary A.l. For any (p,r) G [l,oo] 2 with r hnite, there exists a constant C so that for 
all fi > 0, 

C _1 ||/|| 2 < ^l|e"* A /||rrmi .mrc<^ < C\\f\\ 2 
Besides, for all T > 0, we have 

||e^ A /|| 2 <C||/|| 2 , 

B p /(R d ) B p J(R d ) 

and, if in addition p < 00 then the map T 1— > e^ / is continuous on M.+ with values in 
Blp(R d ), whenever f is in s|~ f (M d ). 



Proof. The first item corresponds to the previous proposition with s = —2/r. Given that 
(e M )j>o is a contracting semigroup over L P (R ), we readily get the second item. The 

■ 2— - 

continuity result is a consequence of the density of smooth functions in B P)T r (R ) if both p 
and r are finite. □ 

To prove Theorem 1.1, we also need the following result: 
Lemma A.l. The operator A defined by 

(A.4) A : f 1— > [t »-> f V 2 e^ (t - T)A / dr} 

is bounded from L r (0, T; L p (R d )) to L r (0, T; L p (R d )) for every T € (0, 00] and 1< p, r < 00 . 
Moreover, there holds 



L l W-Af\\L! r (LP(R d )) - ^ , H/llL5,(LP(R d )) 



and for all T > 0, 

M 1 



f e^ T - t)A fdt 2 2 <C\\f\\ r r (rp(Md)V 
Furthermore, the map T 1— > L e^ *' f dt is continuous on R + with values in B p< 
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Proof. The first part of the statement is Lemma 7.3 of [16]. For establishing the second 
inequality, we just have to notice that 

/ ePV-^fdt 2 2 <C\\Af(T)\\ 2 . 

Using that (e^ )t>o is contracting over L P (R ), and Corollary A.l, we may write 

\\Af(T)\\ 2 ,<C^\\Af\\ Lr(R+ . tLPiRd)y 

By changing / to /1[ot]> one § e ts the desired inequality. The continuity result follows by 
density. □ 

From now on, to simplify the presentation, we agree that A p denotes either the Stokes 
operator on L P (Q) with SI a C 2 bounded domain (see just above Definition 7.1), or the heat 
operator on L P (R ). Lemma A.l extends as follows: 

Lemma A. 2. Let 1 < p, r < oo and Operator A p be defined by 



•Ap '■ J 



t^ [ A p e- {t ~ T)A vf(T)dT ■ 
Jo 

Then for all real number a € [0, 1 — 1/r) there exists a constant C so that for all T > 0, 

II^Aj/IIl^Lp^)) < C\\t a f\\ L r T ( LP ( n) y 

Proof. As regards the heat semigroup in M , this result has been established in [14]. We here 
propose another proof that also works for the Stokes semigroup in bounded domains (and, 
more generally, whenever maximal regularity estimates are available). Let 

(A.5) v(t) d = f e-^ T ^f{r)dT. 



■Jo 
Because 

d t (t a v) + A p (t a v) =t a f + at a ~ 1 v and (t a v)\ t=0 = 0, 
we readily have 

(A.6) \\t a A P v\\ L r T(LPm < C(\\t a f\\ L r T{LP{n)) + a|r-Ml^(Lp(n)))- 

From the definition of v and the fact that (e p )a>o is contracting on L P {Q), we infer that 

K*)||LP(n)< / Wf(r)hnn)dT. 
Jo 

Therefore, 

\\t a -MmL H ^)< [{^f nr)^ with nr) = \\r a f(r)\\ LP(n) , 

that is to say, 

l 



*°-Mt)||LPfm < / (r'r a F(tr')dT'. 
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Hence taking the norm in L r (0,T), and using Minkowski inequality and a + 1/r < 1, 

»1 / r-T 



\t a - l v\\Lr T(L v m < /( r ')~ Q (7 F r (tr')dty dr' 

< I ( T ')- a - l / r ( f T F r (t') dA " dr' 



< c||-^ , IIl i -(o,t)- 

Plugging this inequality in (A. 6) completes the proof. □ 

Remark A.l. Applying the above result in the Stokes case implies that the function v 
defined in (A. 5) and the corresponding gradient term VII satisfy 

(A.7) \\t a (d t V, V 2 V, Vn)|| L r (LP(n)) < C\\t a f\\ L r T(LPm . 

The next lemma provides estimates on the gradient. 
Lemma A. 3. Let 1 < p, q, r < oo with p < q. Let Operator B be defined by 



B:f 



*-)• f Ve- {t - T ^f{ T )dT 
Jo 



(1) If in addition < 1 then for all real numbers a and (3 satisfying 

(A.8) /3 = a + -( ] and ar < 1 

2\p q) 2 

there exists a constant C so that for all T > 0, 

(A-9) \\t $ Bf\\ L r T{Lqm < C\\t a f\\ L r T(LPm . 

(2) // in addition < 1 then for all real numbers a and (3 satisfying 

(A.10) P = a + -( ) h - and ar' < 1 

2 \p q) 2 r 

there exists a constant C so that for all T > 0, 

(A.H) \\t p Bf\\ L¥(Lqm < C\\t a f\\ L r T(LPm . 

(3) More generally, for any s € [r, oo], if < 1 — - + - and (a,/3) satisfy 

(A.12) j3 = a + -{ ) 1 and ar' < 1 

2 \p qj 2 r s 

i/ien i/iere exists a constant C so that for all T > 0, 

(A-13) \\t P Bf\\ L s T{Lqm < C\\t a f\\ L r T(LPm . 

Proof. The limit case /3 = a of the first inequality is a consequence of Lemma A. 2 and 
Sobolev embedding. To treat the case ft < a, the starting point is the following inequality 

(A.14) ||Ve-^/|| L , (n) < Cr s \\f\\ LP{Q) with 5 d ^ f ±(± - i) + \ 

which holds true whenever 1 < p < q < oo. It has been proved in [13] for the Stokes operator 
in bounded domains, and follows from an explicit computation for the heat operator in R . 
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This inequality obviously implies that 



(A.15) tP\\Bf(t)\\ Lq{ a)<CtP f (t-r)- 5 T- a F(r)dT with F(r) d M \\r a f(r)\\ LP(n) . 

Jo 

Therefore, making a change of variables, and using the relationship between a and /3, 

ABf(t)\\ Lq{n) < C f\l-rr 5 {r'r a F{T't)dr'. 
Jo 

Now, taking the L r T norm of both sides and applying Minkowski inequality implies that 

\\t Bf\\ L r T{Lqm <C J\\- T'T 5 {T')- a U F r (r't)dt\ T dr'. 

Then arguing exactly as in the proof of the previous lemma, we get Inequality (A. 9) whenever 
(A. 8) is satisfied. 

In order to establish (A. 11), we start from (A.15) and apply Holder inequality. We obtain 
for all t G [0,T], 



^\\Bf(t)\\ Lq{Q) < Cifi (J\t - r)- r ' s r- ar ' dr^j * \\F\ 



L r (0,T)- 



Then making the usual change of variable and taking advantage of the relationship between 
a and /3, and of the definition of F, we get for all t G [0,T], 



,t 



A\Bf{t)\\L«<si) < c(J {l-rr r ' & {r')- ar 'Y\\t a f\\L^ m - 

It is now clear that we get Inequality (A. 11) under the constraints of (A. 10) and 5r' < 1. 

In order to treat the general case, we have to combine the methods for proving (A. 9) and 
(A. 11). Starting from (A.15), we write 

where </? is a parameter in [0, 1], to be fixed hereafter. 

Applying Holder inequality to the inner integral, we get (with obvious notation) 

I s <C f t 0s ( f (t- T)- S(p rF(r)T- a dr) (I if- t)-^-^)' F(r)r- a dr) dt- 

Applying again Holder inequality, in the last integral only, we thus find that the above r.h.s. 
is bounded by 

s—r 

{t _ T) S(l-myr' T -«r 'dX^dt. 

Then we perform the same change of variables as above to get 

T / r-l \ r / f\ 



r T i 


' /•* » \ r / 


' r* \ r^ 1 / 


' /•* 


1 iP'l 


/ (t-T)- 5 ^F(T)T- a dr) I 


/ F^T)dr) 


/ 


o V 


Jo J V 


Jo J V 


Jo 



( i (1 - t)- s ^ F(tr)T- a drX ( I (l-r)-'( 1 -*')(r) v r- ar 'dr > ) '' 



I s < C\\F\\ s - r ^ T) I II (l-Tr'^~F(lT)r- n dr ] ( / (1 - r )-<HW.H 7 r<- r -,,r (h j (]j 
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where we have used the fact that f3 = a + 5 — 1/r' — 1/s. If we assume that ar' < 1 and that 
5(1 — ip)(-)'r' < 1 then the last integral is bounded. Then applying Minkowski inequality to 
swap the integrals on [0, T] and on [0, 1], we eventually get 

/ < C||F||^ T) (J\l - T )- 5 ^rr- a (J F r (tr) dt) ^rj 8 , 

whence 

\\t^f(t)\\ L s T(Lqm < C||F||^ )T) (J\l - r)- 5 ^r- a --r (^f T J F r (t') dA T dr 

which implies Inequality (A. 13) provided a + 1/r < 1 and 5<ps/r < 1. In order to complete 
the proof, it is only a matter of taking ip = , r so that the conditions 5(1 — <p)(-)'r' < 1 
and 5(ps/r < 1 are equivalent. □ 



Finally we need a lemma involving the following operator 

C:/- 



t-> f e- (t - T)A ?f(T)dT 
Jo 

Lemma A. 4. Let 1 < p, q, r < oo with q > p. 

(1) If in addition - — - < 2 then for all real numbers a and 7 satisfying 

,. x dfl 1\ , , 

(A.16) 7 = a + - -1 and ar'<l 

2\p q) 

there exists a constant C so that for all T > 0, 
(A-17) l|t 7 C/||^( L9( n)) < C\\t a f\\ L r T(LPm . 

(2) If in addition - — - < 2 — - then for all real numbers a and 7 satisfying 

(A.18) 7 = a + -( ) and ar' < 1 

2 \p q) r 

there exists a constant C so that for all T > 0, 
(A.19) ||t 7 C/|| L?? ( L9 (n)) < C\\t a f\\ L r T(LPm . 

(3) More generally, for any s € [r, 00], if < 2 — - + - and (a, 7) satisfy 

(A.20) j = a+-( ) - H and ar' < 1 

2 \p q) r s 

then there exists a constant C so that for all T > 0, 

( A - 21 ) ll i7 C/|Uj,(L9(n)) < C\\t a f\\L r T (LP(n))- 

Proof. The limit case 7 = a of the first inequality is a consequence of Lemma A. 3 and 
Sobolev embedding. To treat the case 7 < a, the starting point is the following inequality 



^ tAvf] \mn)<Ct- & \\f\\ LP(n) with ^ = f ^ (^ - ^ 



(A.22) \\e- tA *f\\Lm)<Ct- 6 \\f\\ Lvm with 5' 



whenever 1 < p < q < 00, which has been proved in [13] for the Stokes operator in bounded 
domains, and follows from an explicit computation for the heat operator in R . 
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This inequality obviously implies that 

(A.23) t~<\\Cf(t)\\ Lq(n) <Ct~< f (t-r)~ 5 T- a F{r)dT with F(r) ^ \\r a f(r)\\ LP{Q) . 

Jo 

Therefore, making a change of variables, and using the relationship between a and 7, 

^\\Bf(t)\\ Lq{a) <C f\l-T')- 5 {r')- a F{T't)dr'. 
Jo 

Now, taking the L r T norm of both sides and applying Minkowski inequality as in the above 

lemmas yields Inequality (A. 17). 

In order to establish (A. 19), we start from (A.23) and apply Holder inequality. We obtain 
for all t £ [0,T], 

* 7 l|B/(t)|| L , (ft) < Cf(j\t-Tr' 5 T- ar ' drY \\F\\ Lr ^ T) . 

Then arguing as in the previous lemma easily leads to (A. 19) under the constraints of (A. 18). 
The general case s > r follows from similar arguments. The details are left to the reader. □ 

We finally want to establish decay estimates for the free solution to heat equation in the 
whole space. 

Lemma A. 5. Assume that no € Bp r (Jl ) with 1 < p, r < 00. The following inequalities hold 
true: 

(1) If s<2 then 

(A. 24) ||t a V e - uo||j:r( R+;L p(] R d)) < C , |ko|| B ^ r ( ]R d) 

(2) If s<l then 

(A. 25) \\t Ve ^o|lLr( R+ . L p( R d)) < C||«o|lB| ir (R d ) 

(3) If s <0 then 

(A. 26) ||t 7 e ^o|lLr( R+;LP ( R d)) < C|l u o|| B | r ( R d) 

Proof. The assumption ensures that V 2 uq € Bp~ 2 (R ). 
yields 

. . 1 S 1 j. A r\ .. . . I S . . 1 A i-y .. .. r\ 

\\t 2 re V U0\\ L r( m+ . L p( R d^ = \\t 2||e V U \\ LP ( R dj\\ L r( R+ .M^ \\V Uollgs-^^dy 

The proof of the other inequalities is totally similar. □ 

Remark A. 2. In this paper, we mainly consider the case where s = —l-\ The corre- 
sponding values of (a, /?, 7) are 

3 d 1 ., d o 1 d 1 ■, d 1 d 1 ., 

a = - - if p > -, 13 = 1-- if p> -, 7 = - - if p>d. 

2 2p r 3 2p r 2 2 2p r 
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